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1»V  ABSTRACT 

Department  of  the  Navy 

Measurements  of  the  spectrum  of  ocean  surface  reradiated  acoustic -signals  generated  by  an 
omnidirectional,  sinusoidaTsource  are  used  to  esti mdfe~rne  ocean  surfoce  wave-height  spectrum. 
Plane-wave  physical  optics  theory  is  applied  to  a  long-crested  sea  surfoce  model  to  verify  Marsh's 
hypothesis  that,  fc^r  low  sed  states,  the  spectrum  of  acoustic  wcves  reradiated  by  the  rough  moving 
sea  surface  consists  of  a  specular  component  at  the  transmitted  frequency  ond  a  sca’.tered  component 
which  is  the  weighted  surface  wave-height  spectrum.  The  approach  employs  theories  by  Parkins, 
Clay  and  Medwin,  and  Pierson.  The  weighting  function  applied  to  the  surface  spectrum  depends 
on  signal  frequency,  experimental  geometry,  ond  mean-square  wave  height.  CW  signals  at  127  Hz 
and  1702  Hz  were  transmitted  from  fixed,  omnidirectional  sources  to  receivers  4.5  Pm  and  41  km 
downrange.  Additionally,  wave  height  was  monitored  continuously  by  resistance  wave  staff.  The 
time  series  records  of  sound  pressure  level  and  ocean  wave  height  were  analyzed  for  spectral 
content.  Agreement  between  theory  and  experiment  is  good  when  the  acoustic  sideband  levels 
are  not  more  than  25  dB  below  the  peak  level  of  the  specular  component.  / 
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ABSTRACT 

An  experimental  and  theoretical  study  demonstrated  that  measurements  of  for* 
ward- scattered  underwater  acoustic  waves  con  be  used  to  estimate  the  ocean  surface 
wave- height  frequency  spectrum,  under  surface  duct  propagation  conditions. 

The  experimental  phase  of  the  investigation  consisted  of  transmitting  127-Hz 
and  1702-Hz  continuous  wave  acoustic  signals  from  fixed,  omnidirectional  sources 
in  the  BiFI  acoustic  range  in  Block  Island  Sound  to  receivers  4.5  km  and  41  km 
downrange.  During  the  acoustic  transmissions,  ocean  wave  height  was  monitored 
continuously  by  means  of  a  resistance  wave  staff.  The  time  series  records  of  sound 
pressure  level  and  ocean  wave  height  were  analyzed  by  digital  computer  to  deter¬ 
mine  spectral  content.  Careful  consideration  was  given  ro  optimizing  signal 
processing  to  achieve  maximum  resolution  and  stability.  The  results  showed  that  the 
surface  reradiated  signal  spectrum  is  spread  in  frequency  and  consists  of  (1)  a 
broadened  specular  component  around  the  transmitted  frequency  and  (2)  sidebands. 
The  sidebands  resemble  the  surface  height  spectrum  over  typical  bandwidths  of 
0. 2  Hz  and  1.0  Hz  for  the  127-Hz  and  1702-Hz  signals,  respectively. 

in  the  theoretical  phase  of  the  investigation,  generalized  physical  optics 
theory  was  used  to  predict  the  spectrum  of  a  sinusoidal  signal  transmitted  by  an 
omnidirectional  source  and  reradiated  by  the  rough,  moving  sea  surface.  This 
generalization,  which  is  based  on  theories  by  Parkins  and  Cloy  and  Medwin,  was 
accompl ishea  by  subdividing  a  large  illuminated  area  on  the  sea  surface  and 
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applying  the  plane  wave  Helmholtz  integral  theory  to  each  subsection.  The  indi¬ 
vidual  autocovariances  were  evaluated,  summed,  and  Fourier  transformed  to  obtain 
the  fdrfield  acoustic  spectrum.  The  employment  of  the  linear,  long-crested  model  of 
random  gravity  waves,  due  to  Pierson,  facilitated  the  integration  of  the  scattering 
equations  for  low  sea  state  conditions.  The  long- crested  waves  are  assumed  to 
propagate  in  a  direction  parallel  to  the  source- receiver  line. 

The  farfield  acoustic  spectrum  is  shown  to  consist  of  a  specular  component  that 
is  free  from  Doppler  at  the  transmitted  frequency  and  a  scattered  component  that 
appears  as  sidebands.  The  sidebands  are  symmetrical  about  the  carrier.  The  distri¬ 
bution  of  the  variance  in  the  sidebands  is  given  by  the  weighted  surface  height 
spectrum.  The  weighting  function  depend:  on  the  signal  frequency,  the  experimental 
geometry,  and  the  mean-square  wave  height.  The  weighting  function  behaves  like 
a  low-pass  filter  that  is  relatively  flat  out  to  some  frequency,  f  j ,  and  then  slopes 
down  to  a  "cutoff"  frequency,  f^  .  Information  about  the  sea  surface  wave-hei^it 
spectrum  can  be  obtained  wily  between  these  two  frequencies.  The  effective  band¬ 
width  of  the  model  at  127  Hz  and  1702  Hz  is  0. 15  Hz  and  0. 5  Hz,  respectively. 

Agreement  between  theory  and  experiment  was  good  when  the  acoustic  side¬ 
bands  were  not  more  than  25  dB  below  the  peak  level  of  the  coherent  component. 
The  bandwidth  of  the  theoretical  weighting  function  was  seen  to  be  less  then  the 
bandwidth  of  the  weighting  function  computed  from  field  data.  This  discrepancy 
can  be  explained  by  noting  that  the  long-crested  sea  surface  model  is  deficient  at 
the  high  frequencies  when  the  sea  surface  is  really  short- crested.  It  appears  that, 
when  the  short- crested  character  of  the  sea  surface  is  included  in  the  scattering 
model  using  methods  due  to  Pierson,  the  results  will  approach  Marsh's  prediction 
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thot  the  first-order  acoustic  sidebands  are  equal  to  the  uniformly  weighted 
surface  height  spectrum  for  unrestricted  geometry. 
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1.0  INTRODUCTION 

Remote  sensing  of  the  environment  offers  many  exciting  possibilities  for 

discovery  in  the  field  of  geophysics.  Aerial  radiometer  flights  to  determine  sea 

surface  temperature  are  now  a  standard  technique  of  physical  oceanography  and 
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marine  meteorology.  Chia  and  Pierson  and  Moore  have  shown  that  it  is  possible 

to  obtain  detailed  information  about  the  ocean  surface  wave  spectrum  and  other 

characteristics  of  the  air-sea  boundary  layer  by  means  of  radar  waves  transmitted 

and  received  by  on  orbiting  space  vehicle.  Underwater  sound  waves  have  been 

used  since  the  1920‘s  to  map  the  ocean  floor  and  to  chart  the  sublayering  of  the 

3 

sediments.  VigoureuxandHersey,  however,  have  noted  that  it  is  unfortunate  that 

few  applications  of  underwater  sound  have  been  made  in  physical  oceanography. 

4 

Clark  and  Yamal  I  have  used  underwater  acoustics  to  study  experimentally  the 

tides,  currents,  and  internal  waves  in  the  Straits  of  Florida.  They  have  suggested 

that  the  underwater  acoustic  range,  consisting  of  a  fixed  source  and  fixed  receivers, 

could  be  made  into  an  effective  remote  probe  to  study  the  marine  environment. 

This  investigation  attempts  to  stimulate  the  interest  of  physical  oceanographers 

in  underwater  acoustics  by  demonstrating,  both  theoretically  and  experimentally, 

that  the  low-frequency  end  of  the  ocean  surface  height  spectrum  can  be  estimated 

from  the  time  history  of  received  sound  pressure  levels.  The  general  problem  of 

using  urxterwater  sound  to  obtain  information  about  the  geometry  and  kinematics 
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of  the  sea  surface  was  probably  first  considered  seriously  by  Eckart  in  1953.  He 


studied  the  cose  of  a  highly  directional,  single- frequency,  acoustic  projector 
transmitting  in  the  direction  of  the  stochastic  sea  surface.  The  transmitter  induces 
secondary  sources  on  the  surface  that  determine  the  nature  of  the  scattered  field  in 
the  water.  The  acoustic  field  at  a  subsurface  receiver  is  then  represented  by 
means  of  the  Helmholtz^  integral,  which  is  evaluated  over  theinsonified  area  on  the 
surface  in  terms  of  the  boundary  values  of  the  pressure  field  and  its  normal  deriv¬ 
ative.  (This  approach  is  often  referred  to  as  the  "physical  optics"  method). 

Eckart  determined  the  scattering  coefficient  for  low  surface  roughness  in  terms  of 
the  directional  surface  wave  spectrum.  He  concluded  that  it  would  be  necessary 
to  vary  the  direction  and  frequency  of  the  acoustic  waves  to  obtain  an  estimate  of 
the  surface  height  spectrum. 

Eckarfs  paper  stimulated  interest  in  the  acoustic  community,  and  c  number  of 
papers  followed  in  which  the  basic  physical  optics  theory  was  modified  and 
extended  in  order  to  make  application  to  more  general  situations  possible.  Proud, 
Beyer,  and  Tamarkin^  used  Eckart's  theory  to  estimate  the  spatial  correlation 

Q 

function  of  a  rough,  reflecting  surface  in  the  laboratory.  Clay  extended  the 

theory  to  include  omnidirectional  sources  by  subdividing  the  insonified  area  on  the 

surface  and  applying  the  Helmholtz  integral  to  each  subdivision  and  summing  the 
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contributions.  He  compared  his  theory  with  sea  data  taken  by  Brown  and  Ricard 

.  10 

and  attempted  to  estimate  the  correlation  distance  on  the  ocean  surface.  Medwin 
showed  that  the  mean-square  wave  height  for  cases  of  small  surface  roughness  can 
be  determined  from  measurements  of  specularly  reflected  acoustic  intensity. 

Other  important  contributions  were  mode  by  Liebermann,  who  introduced 
the  analogy  between  the  reflection  of  sound  waves  at  the  sea  surface  and  the 


reflection  of  plane  waves  by  a  diffraction  grating,  and  by  Marsh,  who  formu¬ 
lated  the  scattering  problem  in  a  way  completely  different  from  that  of  Eckart  and 
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obtained  an  exact  solution  to  the  scattered  field  by  using  Wiener's  theory  of 

generalized  harmonic  analysis.  Marsh's  approach  was  basically  a  randomization 

14 

of  Rayleigh's  solution  to  the  scattering  problem  by  Fourier  series  methods. 

15 

Marsh  and  Mellen  predicted  that  the  first-order  spectrum  of  single¬ 
frequency  acoustic  waves  reflected  from  a  rough,  moving,  isorropic  sea  surface 
should  consist  of  two  distinct  parts,  i.e. ,  to  first  order.  The  first  part,  or  coherent 
component,  is  a  line  at  the  transmitted  frequency,  or  "carrier,"  that  corresponds 
to  a  specular  reflection  from  the  mean  sea  surface.  The  second  part  is  a  continu¬ 
ous  spectrum  produced  by  incoherent  scattering  that  is  symmetric  about  the  carrier 
and  is  a  uniformly  weighted  reproduction  of  the  sea  surface  height  spectrum. 

Parkins^  applied  physical  optics  theory  to  determine  the  spectium  of  the  far- 
field  acoustic  pressure  due  to  an  incident,  highly  directional,  single- frequency 
signal  reradiated  from  a  rough,  moving  surface  described  by  the  Neumann-fierson 
directional  wave  spectrum.  The  received  signal  spectrum  is  shown  to  consist  of  a 
line  component  at  the  carrier  frequency  (the  coherent  component)  and  two  delta 
functions  positioned  at  frequencies  equidistant  from  the  carrier.  The  location  of 
the  delta  functions  depends  on  the  angles  of  incidence  and  reception  relative  to 
the  wind  direction  and  on  the  frequency  of  the  transmitted  signal.  The  magnitude 
of  the  delta  functions  is  determined  by  the  surface  wave-height  spectrum.  This 
result  does  not  necessarily  contradict  the  predictions  of  Marsh  and  Mellen,  because 

Perkins  considers  a  highly  directional  source,  whereas  they  allow  for  a  wide  range 
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of  angles  of  incidence  and  reflection.  Gulin  and  Roderick  and  Cron  adopted 
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physical  optics  theory  to  explain  scattering  from  moving,  sinusoidal  surfaces. 

19 

Nuttail  and  Cron  used  physical  optics  theory  to  derive  the  spectrum  for 
highly  directional  signals  rjradiated  from  narrow- band,  time- varying  surfaces 

described  by  particular  correlation  functions. 

20 

Clay  and  Medwin  extended  Clay's  earlier  theory  to  determine  the  covari¬ 
ance  of  signals  reradiated  from  a  slightly  rough,  moving  water  surface  described  by 
a  traveling  correlation  function. 

Experimental  verification  of  these  theories  has  been  extensive  in  the  labora¬ 
tory,  but  relatively  few  comparisons  hav  been  made  with  data  taken  at  sea, 

especially  for  the  case  of  forward  scattering. 

21 

Mackenzie  reported  on  a  series  of  experiments  in  which  pure  tones  of 
frequencies  350  Hz  to  2400  Hz  were  transm. . .-ed  across  a  shallow  water  area  and 
detected  at  ranges  of  up  to  27  km.  The  received  signal  spectrum  was  spread  in 
frequency  in  such  a  manner  that  the  relative  power  (Pr)  was  related  to  the  half- 
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width  of  the  spectrum  J  f  -  fg  by  Pf  =  b^fg  f  -  fg  ,  where  b  is  a  constant, 

fg  is  the  transmitted  frequency,  and  f  is  the  measured  frequency  at  the  receiver. 
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Scrimger  measured  the  phase  and  amplitude  fluctuations  produced  by  waves 
generated  by  a  boat  running  between  a  source  and  receiver  separated  by  a  distance 
of  100  m.  He  found  that,  for  slightly  rough  surfaces,  the  frequency  spectrum  of 
the  acoustic  fluctuations  was  the  same  as  the  frequency  spectrum  of  the  water 
surface.  For  rough  water  surfaces,  the  signal  spectrum  was  broader  and  exhibited 


second  and  third  harmonics. 


Gulin  and  Malyeshev  measured  the  spatial  correlations  of  the  amplitude  and 


phase  fluctuations  of  signals  reflected  from  the  surface  of  a  fresh  water  reservoir. 


They  found  that  “the  signal  amplitude  spatial  correlation  coefficients  in  the  case 

of  small  values  for  the  Rayleigh  parameter  usually  take  the  form  of  damped- 

cscil lotion  functions. . .  .The  oscillation  period  is  related  to  the  average  wave** 

length  of  these  water  waves  on  the  sea  surface,  in  the  case  of. . (large  Rayleigh 

parameter)  a  more  abrupt  drop  in  the  correlation  coefficient  is  observed  and  the 

quasiperiodic  behavior  disappears.  “ 
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Clark  and  Yamal  I  discussed  an  experiment  performed  at  Fowey  Rocks  Light¬ 
house,  Miami,  in  which  an  acoustic  continuous  wave  (CW)  signal  was  reflected 
from  the  sea  surface  and  received  by  a  hydrophone  137  m  away  in  18  m  of  water 
(Figure  1).  The  resemblance  between  the  wave-height  power  spectral  density  (PSD) 


Figure  1.  “Single-reflection"  surface  wave  experiment 
(From  Clark  and  Yarnall.^) 


and  the  phase  and  amplitude  PSD  is  evident. 

24 

Ellinthorpe  transmitted  low-frequency  signals  from  a  fixed  source  to  a  fixed 
receiver  located  in  the  Sargasso  Sea.  He  observed  that  "the  spectra  (Figure  2)  all 
show  deep  nulls  on  either  side  of  the  transmitted  frequency  and  a  pair  of  relatively 
broad  sidelobes  centered  at  approximately  0. 1  Hz  away  from  the  central  peak.  . . . 
These  results  we  consistent  in  several  respects  with  what  would  be  expected  if  the 
underlying  physical  process  were  a  phase  modulation  imposed  on  the  acoustic  wave 
by  reflection  from  the  moving  ocean  surface." 
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Figure  2.  Received  spectrum  of  a  surface -reradiated, 
single- frequency  acoustic  signal 
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Roderick  and  Cron  measured  surface  reflected  sound  in  the  deep  ocean  near 
Bermuda  and  compared  the  spectral  content  of  the  received  signals  with  the  surface 
wave-height  amplitude  spectrum  measured  25  nmi  away  at  Argus  Island.  The 
results  shewed  that  new  frequency  components  in  tf  e  acoustic  signal  were  gener¬ 
ated  by  the  sea  surface  reflection  process  and  that  these  components  weredisplaced 


from  the  carrier  at  multiples  of  the  peak  frequency  in  the  suiface  height  spectrum. 

The  experiments  cited  above  are  of  considerable  value  because  they  show  that 

single- frequency  acoustic  signals  reradiated  from  a  rough,  moving  water  surface 

are  spread  in  frequency  and  that  this  spread  is  closely  related  to  the  properties  of 

the  surface.  It  therefore  seems  reasonable  that  the  ocean  wave-height  frequency 

spectrum  can  be  inferred  from  the  acoustic  measurements.  The  purpose  of  this 

thesis  is  to  show  by  theory  and  experiments  that  the  ocean  wave-height  spectrum 

can  be  determined  explicitly  from  the  acoustic  forward -scattered  spectrum.  The 

theoretical  method  used  is  a  combination  of  the  theories  of  Parkins^  and  Clay  and 
20 

Medwin.  The  theoretical  predictions  are  then  compared  with  spectra  estimated 
from  field  measurements  of  sound  pressure  level  and  ocean  wave  height  taken  in 
Block  Island  Sound.  Hopefully,  methods  similar  to  those  developed  here  could  be 
applied  to  obtain  information  on  ocean  currents,  internal  waves,  and  turbulence. 
Before  considering  a  particular  theoretical  formulation  of  the  scattering  problem, 
however,  a  brief  review  of  pertinent  underwater  sound  transmission  characteristics 
is  presented. 


2.0  CHARACTERISTICS  OF  UNDERWATER  ACOUSTIC 
TRANSMISSIONS  IN  THE  SEA 


The  science  of  ocean  acoustics  deals  with  the  effects  of  the  ocean  environment 
on  underwater  sound  propagation.  The  sources  of  sound  in  the  sea  are  either 
natural  or  man-made.  Examples  of  natural  sounds  are  the  cracking  of  shrimp,  the 
communications  of  porpoises  and  whales,  the  action  of  the  wind  on  the  waves,  and 
seismic  disturbances.  Examples  of  man-made  sounds  are  the  explosions  used  in 
seismic  studies,  signals  from  underwater  communication  systems  and  sound  detection 
and  ranging  systems  (sonar),  and  cavitation  from  the  propellers  of  ships.  In  this 
study,  only  sinusoidal  signals  between  100  and  2000  Hz  will  be  considered. 


2. 1  THE  VELOCITY  OF  SOUND  IN  THE  OCEAN 

AND  ITS  EFFECTS  ON  ACOUSTIC  TRANSMISSION 

The  propagation  of  sound  waves  is  described  to  a  high  order  of  approximation 
by  the  scalar  wave  equation: 


V2<X>  = 


1  d2(f> 

T  T 
c  at 


0) 


where 

<t>  =  acoustic  pressure 

C  =  speed  of  sound  expressed  in  Cartesian  coordinates  x,  y,  and  z. 

When  the  variations  in  C(x,  y,  z)  are  known,  the  acoustic  field  intensity  can  be 
determined,  in  principle,  by  means  of  Eq.  (1).  In  most  studies  of  sound 
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propagation  over  relatively  short  ranges,  only  vertical  variations  in  C(x,  y,  z) 

have  been  considered  (i.  e. ,  C(x,  y,  z)  «  C(z) )  .  The  speed  of  sound  in  ?ea  water 

is  a  function  of  the  density  and  the  isentropic  bulk  modulus;  it  ran  be  determined 

from  routinely  measured  oceanographic  quantities,  i.e. ,  temperature,  salinity,  end 

depth.  Several  theoretical  and  empirical  relationships  have  been  derived  that  give 

C  as  a  function  of  these  variables.  The  most  generally  accepted  equation  was 
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obtained  by  Wilson  as  a  polynomial  best  fit  to  experimental  data  measured  in  the 
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laboratory.  This  equation  is  given  by  Kinsler  and  Frey  as  follows: 

C  =  H+9+4.6T  -  0.055T2  +0t0003T3+  (1 . 39  —  0.01 2T)  (S  -35)+0.017h,  (2) 

o 

where 

C  *  speed  of  sound  in  meters  per  second 
T  =  temperature  in  degrees  Celsius 
Sq  =  salinity  in  parts  per  thousand 
h  =  depth  in  meters. 

From  typical  oceanographic  data,  this  equation  shows  that  a  1°  C  rise  in  temper¬ 
ature  increases  the  sound  speed  by  about  3  m/sec,  a  increase  in  salinity 
increases  the  sound  speed  by  1.  3  m/sec,  and  10  m  increase  in  depth  increases  the 
sound  speed  by  0. 17  m/sec. 

Since  the  ocean  is,  in  general,  stratified  rather  than  homogeneous,  the  speed 
of  sound  is  a  function  of  depth.  This  results  in  acoustic  energy  being  refracted 
toward  regions  of  minimum  velocity  as  a  consequence  of  Fermat's  principle,  which 
reduces  to  Snell's  Law  for  C(x,  y,  z)«C( z)  .  The  speed  of  sound  increases  with 
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increasing  <,<spth  in  the  near-surface  layer  of  the  Arctic  and  Antarctic  oceans,  in 

the  mixed  layer  of  moderate  latitudes,  and  in  the  coastal  waters  of  temperate  zones 

in  winter.  Within  the  near  surface  layer,  propagating  sound  waves  are  refracted 

upward  toward  the  region  of  lower  sound  speed  near  the  sea  surface.  At  the  sea 

surface,  the  large  acoustic  impedance  mismatch  between  water  ur>d  air  causes  a 

nearly  complete  reflection  of  acoustic  waves.  The  reflected  waves  are  bent  by 

refraction  toward  the  sea  surface,  and  the  process  is  repeated  again  and  again, 

thus  fulfilling  the  condition  for  a  surface  channel  or  duct.  The  duct  affords  a 

means  of  iong-renge  acoustic  signal  transmission  with  small  associated  energy  loss. 

27 

The  important  features  of  these  ducts  are  discussed  by  Schulkin. 

Acoustic  signals  reflected  one  or  more  times  from  the  ocean  surface  experience 
energy  loss  as  well  os  amplitude  and  phase  variations.  This  total  energy  loss  is  the 
result  of  the  divergence  of  the  acoustic  waves  from  the  source  (spreading  loss),  the 
reflection  loss  at  the  sea  surface,  and  the  absorption  within  the  water  column  due 
to  viscosity,  thermal  conduction,  and  relaxation.  Interference,  as  well  as  dif¬ 
fractive  leakage  out  of  the  duct,  results  in  additional  attenuation.  The  amplitude 
and  phase  fluctuations  can  be  caused  by  the  rough,  moving  surface  of  the  sec, 
variable  currents,  internal  waves,  end  turbulence.  In  this  study,  only  the  fluctu¬ 
ations  caused  by  (he  moving  sea  surface  will  be  considered.  Signal  fluctuations 
induced  by  the  other  dynamic  processes  in  the  ocean  medium  are  assumed  to  be 
either  small  in  comparison  to  the  surface  wave  effects  or  of  a  different  temporal 
character,  both  of  which  would  result  in  spectra!  signatures  readily  distinguishable 
from  the  sea  surface  effects,  Clark  and  Yarnall^  and  Huff^  have  demonstrated 
that  the  various  effects  of  the  processes  ■  ■>  the  ocean  can  be  distinguished  fronieach 


other  by  time-series  analysis. 
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2.2  ACOUSTIC  PROPAGATION  MODELS 

Ray  theory,  normal  mode  theory,  and  empirical  methods  are  generally  used  for 
the  purpose  of  constructing  models  to  predict  acoustic  field  intensity.  Ray  theory 
is  analogous  to  ray  optics  and,  as  such,  is  intuitively  appealing.  Rays  are  used  to 
show  the  direction  of  propagation  of  acoustic  energy.  If  the  change  in  slope  of  the 
sound  velocity  profile  is  not  too  large  over  a  distance  of  1  wavelength,  the  ray 
solution  is  usually  good. 

In  the  normal  mode  theory,  the  acoustic  field  is  expressed  as  the  sum  of  the 
characteristic  functions  that  satisfy  the  wave  equation  and  the  boundary  conditions. 
The  formulas  generated  are  usually  quite  complicated  and  must  be  evaluated 
numerically  on  a  digital  computer.  The  solution  may  be  difficult  to  interpret  if 
more  than  a  few  modes  are  present  in  the  channel.  However,  in  shallow  water  and 
at  long  range,  the  normal  mode  theory  is  generally  superior  to  ray  theory. 

The  empirical  method  utilizes  large  numbers  of  measurements  of  transmission 

loss  versus  factors  pertaining  to  the  test  geometry  and  the  ocean  medium.  Some 

features  of  the  models  obtained  by  empirical  methods,  particularly  the  Acoustic 

Meteorological  and  Oceanographic  Survey  (AMOS)  model,  are  discussed  by 
27 

Schulkin. 

The  best  choice  for  cm  acoustic  model  is  highly  dependent  on  local  conditions, 
as  well  as  on  what  the  investigator  hopes  to  learn  from  the  model.  In  this  study, 
the  emphasis  is  on  the  reflection  of  acoustic  waves  at  the  sea  surface  without 
regard  to  any  particular  propagation  model,  since  the  purpose  is  to  show  that  the 


IZ 

reflected  signal  spectrum  contains  information  about  the  surface  height  spectrum. 
However,  results  from  both  ray  theory  and  normal  mode  theory  will  be  referred  to 
•n  the  interpretation  of  experimental  data,  as  required. 


2.3  REFLECTIONS  OF  ACOUSTIC  WAVES 
FROM  THE  SFA  SURFACE 

Let  us  assume  a  unit  amplitude  plane  wave  in  the  water: 
$  -  exp  i(aj  x+  b^y  -  c^  z  -ut)  , 

where 

,  b,,  Cj  =  direction  cosines 

x,  y,  z  =  Cartesian  coordinates 
u  =  radian  frequency  =  2trf 
I  =  time  . 


(3) 


If  this  plane  wave  is  incident  on  a  flat  .lir-sea  interface  z  =  0,  then  the 
following  boundary  conditions  must  be  satisfied: 

1.  The  acoustic  pressure  must  be  continuous  across  the  boundary. 

2.  The  \  "ticai  particle  velocity  must  be  continuous  across  the  boundary. 
These  conditions  are  expressed  as 


(5) 
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where 

A,  W  refer  to  air  and  water,  respectively 
p  refers  to  density. 

When  these  boundary  conditions  are  applied  to  the  incident  wave,  the  angle  of 
incidence  is  equal  to  the  angle  of  reflection,  and  the  angle  of  the  refracted  ray 
can  be  determined  from  Snell's  Law.  In  addition,  the  frequencies  of  the  reflected 
and  refracted  waves  are  the  same  as  the  frequency  of  the  incident  wave.  Because 
of  the  large  acoustic  impedance  mismatch  between  water  and  air  >>  (ffc)^)/ 

almost  all  the  acoustic  energy  is  reflected  back  into  the  water.  Thus,  to  a  high 
level  of  approximation  the  reflection  coefficient  is  -  1,  which  means  that  the 
reflected  wave  has  undergone  a  phase  change  of  ir  radians  relative  to  the  incident 
wave.  A  surface  tor  which  (pQ^  >  >  (pC)^  is  called  a  "pressure  release" 
surface.  Thus,  the  sum  of  the  incident  and  reflected  pressures  is  rero: 

^incident +  ^reflected  ®  ® 

Hence,  a  perfectly  calm  sea  surface  is  an  almost  perfect  reflector  of  sound. 

The  problem  of  determining  the  reflected  field  for  a  rough  sen  surface  is  much 
mere  complicated  than  it  is  for  a  flat  sea  surface.  In  general,  an  exact  solution  is 
not  possible.  Physically,  acoustic  energy  is  scattered  in  directions  other  than  the 
specular  direction.  When  the  energy  scattered  by  the  rough,  moving  surface  is 
added  to  the  specularly  reflected  component  of  the  field  (the  coherent  component), 
amplitude  and  phase  fluctuations  are  introduced  into  the  secondary  acoustic  field. 
This  hes  the  effect  of  broadening  the  spectrum  of  the  received  signal,  which  is 
called  "Doppler  broadening"  or"surface  wave  modulation." 
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One  of  the  first  considerations  in  any  study  of  scattering  from  rough  surfaces 

is  to  obtain  an  adequate  definition  of  a  rough  surface.  Perhaps  the  most  common 

definition  of  surface  roughness  is  the  Rayleigh  "roughness  criterion,"  which  states 

that  a  surface  is  rough  if  the  height  of  the  roughness  elements  h  >  — ^ _  , 

8  sin  X 

where  X  is  the  grazing  angle  and  X  is  the  wavelength  of  the  incident  radiation. 

The  Rayleigh  criterion  for  incident  rays  on  a  rough  surface  will  now  be  derived. 

From  Figure  3,  it  is  seen  that 

a  =  h  /sin  X 
r 

b  =  a  cos  0  =  a  cos  (*  -  2X)  .  (7) 

The  path  difference  of  rays  reflected  at  Q  and  at  A ,  Ar  =  a  +  b  ,  is 

h  h 

Ar  =  — —  +  — —  cos  (w  -  2  X) 
sinX  sinX 

or 

Ar  =  A.  [l  +  2sin2X-l]  =2h  sinX  .  (8) 

•  -v  *  jr 

sin  A 


Figure  3.  Derivation  of  the  Rayleigh  Criterion 
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Hence,  when  the  phase  difference  Ar^  is  near  zero,  the  rays  will  be  almost 

in  phase;  however,  if  it  is  near  ir  ,  there  will  be  maximum  destructive  interfer¬ 
ence.  Arbitrarily  choosing  a  value  for  the  phase  difference  that  is  halfway 

between  rough  (*)  and  smooth  (0),  i.e. ,  ir/2  ,  yields  the  Rayleigh  criterion. 

29 

Contemporary  authors  have  suggested  other  criteria,  but  Beckmann  states  that, 
rather  than  develop  a  precise  definition  of  roughness  in  this  way,  it  makes  more 
sense  to  observe  that  a  surface  will  be  smooth  only  if  _L~»  0  or  X -*  0  .  The 


X 

meaning  of  the  Rayleigh  criterion  for  underwater  acoustics  is  that,  as  the  grazing 
angle  decreases  and  the  wavelength  increases,  the  sea  surface  behaves  more  and 
more  i  ike  a  smooth  surface. 

Another  basic  problem  is  to  determine  how  the  scattered  intensity  at  the 
receiver  varies  as  a  function  of  position  on  the  rough  surface.  In  the  case  of 

5 

Eckart's  physical  optics  theory,  a  highly  directional  transmitter  insonifies  a 

relatively  small  region  on  the  surface  around  the  specular  point.  However,  in  the 

case  of  an  omnidirectional  source,  a  very  large  area  on  the  surface  is  insonified 

and  regions  far  from  the  specular  point  may  contribute  significantly  to  the  scattered 
29 

field.  Beckmann  has  examined  this  problem  in  the  case  of  electromagnetic  waves. 
He  regarded  the  reflecting  surface  as  an  aperture  that  diffracts  incident  waves 
propagating  from  a  point  source  to  a  point  receiver.  The  diffraction  pattern 
observed  at  the  receiver  can  be  obtained  by  constructing  Fresnel  zones  on  the 
reflecting  surface.  (These  are  zones  where  the  phase  difference  of  the  radiation 
arriving  at  the  receiver  from  a  given  zone  is  one-half  a  wavelength  different  from 
the  radiation  arriving  from  the  adjacent  zones.)  In  classical  optics,  where  one  is 
concerned  with  the  diffraction  of  light  from  an  aperture,  it  is  shown  that  waves 


16 


from  each  zone  Interfere  destructively  with  those  *rom  the  two  adjacent  zones. 

Consequently,  the  wave  amplitudes  from  ail  zones,  except  the  first  and  the  last, 

nearly  cancel  each  other.  The  total  amplitude  is  then,  to  a  high  order  of 

approximation,  the  average  contribution  from  the  first  and  last  zones.  However, 

since  the  variation  of  the  obliuuity  factor  reduces  the  contrition  of  the  last  zone 

to  almost  zero,  the  result  is  half  the  value  of  the  first  zone.  (The  obliquity-  factor 

determines  how  the  amplitudes  of  the  secondary  waves  vary  with  direction.) 

29 

Beckmann  has  pointed  out  that  caution  must  he  used  in  cpplying  the  results  of 
classical  optics  directly  to  radio  propagation  problems.  He  has,  however,  derived 
formulas  to  describe  the  diffraction  by  the  earth  of  obliquely  incident  radio  waves 
and  has  reached  the  following  conclusions: 

1.  The  first  Fresnel  zone  is  the  most  important  region  on  the 
scattering  surface  in  contributing  to  the  total  field  received  at  a  given 
point. 

2.  Other  important  contributions  come  from  areas  on  the  surface 
above  the  transmitter  and  receiver.  These  areas  are  not  included  in 
any  Fresnel  zc.ne. 

3.  The  Fresnel  zones  on  the  surface  are  ellipses  and,  for  small 
grazing  angles,  are  very  long  and  narrow. 

There  is  every  reason  to  believe  that  these  conclusions  are  valid  for  under- 

22 

water  acoustics.  Scrimger  found  that  water  waves  immediately  above  the 
transmitter  and  receiver  cause  amplitude  and  phase  fluctuations  in  the  received 
signal.  Clay  and  Medwin  and  Clay  found  that  the  areas  far  from  the  specular 
point  must  be  considered  to  correctly  predict  the  scattered  field  intensity. 
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An  important  conclusion  that  has  emerged  from  many  of  the  theoretical  and 
experimental  investigations  on  scattering  from  rough  surfaces  is  that  periodic, 
rough  surfaces  tend  to  behave  as  diffraction  gratings  and  thus  scatter  sound  into 
discrete  directions  corresponding  to  the  orders  of  the  grating  (Bragg  scattering). 

This  means  that  the  directions  of  maximum  intensity  of  the  field  scattered  by  a 
periodic  surface  are  given  by 

sin  ^  =  sin  8  +  p  —  (p  =  0 ,  +  1 ,  +  2 ,  .  .  .  )  ,  (9) 

^  A 

where 

\lr  ~  angle  of  scatter 
P 

0  =  angle  of  incidence 
\  -  acoustic  wavelength 
A  =  surface  wavelength 
p  =  integer  mode  index  . 

29 

This  aspect  of  scattering  has  been  treated  in  detail  by  Beckmann,  who  presents 
scattering  diagrams  as  a  function  of  angle  of  incidence.  In  reality,  cf  course, 
other  directions  of  scattering  are  possible,  but  the  maxima  of  the  scattered  field 
are  g:ven  by  the  Bragg  relation.  Beckmann  has  handled  this  situation  by  removing 
the  constraint  of  the  integer  mode  index  p  . 

Liebermann 1  ^  asserted  that  a  random  surface,  which  is  presumably  composed 
of  elementary  sinusoidal  waves,  behaves  like  a  diffraction  grating  in  the  sense  that 
single- frequency  radiation  will  be  preferentially  scattered  according  to  the  wave¬ 
lengths  present  in  the  surface.  He  demonstrated  experimentally  that,  for  sinusoidal 
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radiation,  a  very  narrow  band  of  frequencies  in  the  rough  surface  produces  the 

scattering.  Therefore,  "a  spectrum  analysis  of  surface  roughness  can  be  obtained 

by  slowly  varying  the  frequency  of  the  incident  radiation  and  observing  the  magni- 

11  12  15  31 

tude  of  the  scattered  radiation!1  Marsh  et  al.,  '  '  in  a  series  of  papers, 

further  developed  this  concept  and  referred  to  an  "equivalent  surface  wave 
number1 1 : 


K  =  k(cosX.  +  cosX  )  , 
in  re 


(10) 


where 


k  =  radiation  wave  number 


X}n  =  incident  grazing  angle 
Xre  =  reflected  grazing  angle  . 
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Wright  employed  this  concept  in  connection  with  the  scattering  of  radar  waves 

from  the  sea  surface.  Recent  evidence  for  the  validity  of  Bragg  scattering  has 

33 

come  from  model  tank  experiments  by  Kingsbury. 

The  concept  of  Bragg  scattering  can  also  be  applied  to  explain  the  Doppier 
broadening  observed  when  a  signal  is  reflected  from  a  rough,  moving  surface.  The 
signal  is  scattered  into  discrete  orders  by  elements  of  the  rough  surface.  Since 
these  elements  are  in  motion,  the  frequencies  of  the  scattered  orders  are  shifted 
away  from  the  carrier.  Gul  in  ^  determined  that  the  frequency  shift  is 


Af  =  p-  , 

A 


(11) 
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wh^re 

v  =  phase  speed  of  the  surface  irregularities 
A  =  wavelength  of  the  surface 
p  =  the  scattered  order  . 

Hence,  there  is  no  Doppler  shift  in  the  specular  direction  where  order  p  =  0  . 
33 

Kingsbury  has  experimentally  verified  this  relationship  for  water  surfaces 
characterized  by  a  sharply  peaked  spectrum. 
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3.0  THEORY 


The  purpose  of  this  chapter  is  to  provide  a  theoretical  basis  for  estimating  the 

ocean  wave-height  spectrum  by  means  of  the  acoustic  forward-scattered  spectrum. 

There  are  many  possible  theoretical  approaches  to  the  general  problem  of  determin- 

34 

ing  the  field  scattered  by  the  rough  sea  surface.  Fortuin  has  recently  provided  a 
fairly  complete  summary  of  the  better  known  methods.  The  majority  of  these 
methods  lead  to  solutions  that  are  not  exact  and  involve  a  large  number  of  simpli¬ 
fying  assumptions. 

Perhaps  the  most  nearly  complete  theory  of  scattering  of  acoustic  waves  in  the 

ocean  is  that  of  Marsh  et  al.  ^  Without  considering  the  time  variation . 

of  the  ocean  surface.  Marsh  demonstrated  by  means  of  a  randomized  Rayleigh -type 

solution  that,  for  small  surface  roughness,  the  received  acoustic  intensity  consists 

of  a  specular  reflection  term  and  a  nonspecular  scattered  term.  Rayleigh-based 

methods  are  frequently  referred  to  as  "wove  expansion  methods. "  The  wave 

expansion  method  is  used  to  derive  expressions  for  the  magnitude  of  the  specular 

and  nonspecular  terms.  These  terms  are  used  in  the  prediction  of  the  acoustic 

surface  loss  per  bounce  in  cases  of  surface  duct  propagation. 

15  31 

Marsh  etal.,  '  through  keen  physical  insight,  have  heuristically  extended 
their  mathematical  model  to  predict  the  effects  of  acoustic  reflections  from  the 
time- varying  surface  of  the  ocean.  They  assert  that  the  first-order,  near-specular 
scattered  spectrum  in  the  frequency  domain  will  be  equal  to  the  uniformly  weighted  sea 
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surface  spectrum  (for  low  surface  roughness).  This  conclusion  was  reached  by 
considering  the  constraints  placed  on  acoustic  propagation  via  surfaced  reflected 
paths,  i.  e. ,  Bragg  scattering,  Doppler  shifts,  the  dispersion  relation  for  deep 
water  gravity  waves,  and  the  requirements  for  the  conservation  of  energy. 

Results  from  simulations  of  scattering  problems  on  the  computer  have  supported  the 
hypothesis  that  the  weighting  is,  in  fact,  uniform  when  ail  possible  angles  of 
scatter  are  considered. 

In  the  present  study.  Marsh's  predictions  of  the  forward-scattered  spectrum  in 

the  frequency  domain  will  be  examined  quantitatively  by  means  of  c  simpler  theory; 

5  38 

i.e.,  the"  physical  optics"  theory  developed  by  Eckart  in  1953.  (Tolstoyand  Clay 

i  / 

develop  this  theory  from  basic  principles.)  Parkins  has  used  physical  optics 
theory  to  determine  the  spectrum  of  a  sinusoidal  acoustic  signal  reradiated  from  a 
rough,  moving  sea  surface  characterized  by  the  Neumann-Pierson  directional  wave 
spectrum.  He  found  that  the  received  acoustic  signal  spectrum  consists  of  three 
delta  functions.  This  result  would  follow  for  any  directional  wave  spectrum  as 
long  as  the  dispersion  relation  holds. 

Parkins'  result  cannot  be  compared  with  Marsh's  predictions,  because  Parkins 
used  a  highly  directional  source  and  receiver  and,  thus,  obtained  a  very  restric¬ 
tive  relationship  between  the  ocean  wave- height  spectrum  and  the  acoustic 
frequency  and  geometry  of  the  experiment,  as  did  Eckart.  However,  it  is  shown 
here  that,  if  Parkins'  theo>y  is  generalized  to  allow  for  wide-beam  sources,  the 
reradiated  acoustic  spectrum  is  equal  to  the  weighted  ocean  surface  wave-height 
spectrum  and  that  the  weighting  depends  upon  the  source  frequency  and  the  experi¬ 
mental  geometry.  Parkins'  theory  is  generalized  for  the  long-crested  sea  surface 
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(Pierson  )  by 
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using  methods  developed  by  Cloy  and  Medwin  and  Nuttoll 

and  Cron.  "  When  the  mode*  is  extended  to  the  short- crested  sea  surface 
39  40. 

(Pierson  '  },  the  weighting  will  probably  be  more  nearly  uniform  and,  hence, 

will  approach  Marsh's  prediction. 

3.1  SYNOPSIS  OF  THE  PHYSICAL 
OPTICS  METHOD 

The  acoustical  pressure  due  to  a  sinusoicjfti,  narrow-beam  signal  reradiated 

from  a  rough,  moving  sea  surface  is  derived  by*|jsing  the  Helmholtz^  integral  and 

the  Kirchoff^  approximation.  This  solution  cjfg  be  generalized  to  describe  an 

8 

omnidirectional  source  by  using  Clay's  method  of  subdividing  the  total  insonified 

38 

area  on  the  surface  and  applying  Helmholtz- Kirchoff  theory  to  each  subarea. 

The  total  field  can  then  be  obtained  by  summing  the  contributions  from  each  sub- 

area,  but  the  coherent  and  incoherent  components  of  the  field  must  be  summed 

separately.  It  is  assumed  that  the  surface  is  rough  in  only  one  dimension  (long- 

crested),  so  that  only  variations  parallel  to  the  source- receiver  line  are  considered. 

Since  the  received  signal  spectruf  is  the  desired  result  of  this  analysis,  the 

autocovariance  function  at  a  single  receiver  due  to  any  two  subareas  (subsections) 

m  and  n  is  formed  from  the  expressions  for  the  received  pressure  by  using  methods 

20 

developed  by  Cloy  and  Medwin.  The  ensemble  average  is  taken  by  using  the 
method  of  characteristic  functions.  ^  When  the  distribution  of  surface  heights  is 
Gaussian,  the  ensemble  average  is  an  exponential  expression  containing  the 
space- time  correlation  function  of  the  sea  surface.  This  exponential  expression  is 
expanded  and  only  the  first  two  terms  are  retained  in  order  to  examine  in  greater 
detail  the  dependence  of  the  acoustic  field  on  the  sea  surface  correlation  function. 


Tit  is  procedure  limits  the  subsequent  analysis  to  cases  of  small  surface  roughness  or 
long  acoustic  wavelengths. 

The  expansion  of  the  exponent  also  separates  the  acoustic  field  contributions 

20 

into  specular  (coherent)  and  scattered  components  and  facilitates  the  eventual 
summation  over  all  subsections.  The  coherent  component  can  be  integrated  and  is 
shown  to  depend  on  the  frequency,  the  geometry  of  the  experiment,  and  on  the 
mean-square  wave  height.  The  scattered  component  is  evaluated  by  using  the 
long- crested,  linear,  Gaussian,  Eulerian  mode!  (Pierson^' ^),  which  expresses 
the  space-time  correlation  function  of  the  sea  surface  in  terms  of  the  ocean  surface 
height  spectrum.  The  resulting  expression  for  the  acoustic  covariance  can  then  be 
integrated  and  Fourier  transformed  by  using  methods  developed  by  Parkins^  to 
yield  the  received  spectrum.  When  the  spectral  contributions  resulting  from  all 
subsections  are  combined,  the  scattered  acoustic  spectrum  is  expressed  as  the 
product  of  the  ocean  wave- height  spectrum  and  a  weighting  factor  ^iat  depends  on 
the  transmuted  frequency,  the  mean-square  wave  height,  and  the  source-receiver 
geometry.  It  is  thus  possible,  in  principle,  to  apply  this  weighting  factor  to  the 
acoustic  scattered  spectrum  to  recover  the  ocean  wave -height  spectrum  if  the 
ocean  surface  is  not  too  rough  and  is  nearly  long- crested. 

3.2  DERIVATION  OF  THE  ACOUSTIC 
surface  RERADiATED  SPECTRUM 

In  the  subsequent  derivation,  a  number  of  approximations  and  assumptions  are 
made.  For  convenience,  they  are  listed  below  as  follows: 

1.  The  acoustic  signal  is  received  in  the  farfield  (Fraunhofer  zone) 


of  hie  source. 


2.  The  plane  wave  approximation  is  valid  for  small  sections  of  the 


insonified  area. 

3.  The  incident  signal  is  a  pure  tone  (sinusoid). 

4.  The  ocean  surface  is  a  pressure- release  surface. 

5.  The  rough  ocean  surface  is  long- crested,  has  a  Gaussian  distri- 
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bution  of  heights,  and  has  a  mean-square  height  of  a  . 

6.  The  surface  is  homogeneous  and  stationary  over  the  time  of  the 
measurement;. 

7.  The  Kirchoff  approximation  is  made  locally,  i.e. ,  the  derivative 
of  the  secondary  field  normal  to  the  surface  is  given  by 


i 


where 

R  =  reflectivity  coefficient,  which  is  -1  for  a  pressure-release  surface 
=  unit  normal  to  the  rough  surface  (positive  down) 

<t>.  =  incident  acoustic  field. 

8.  The  surface  is  not  too  rough,  and  the  reradiated  field  can  be 
separated  into  a  coherent  (specular)  component  and  a  scattered  component. 

9.  The  ocean  medium  is  isovelocity  for  the  depths  considered. 

A  highly  directional  acoustic  transmitter  will  now  be  assumed  in  order  to 

obtain  an  exp.ession  for  the  surface  reradiated  ocoustic  pressure,  os  was  done  by 
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Tolstoy  and  Clay,  who  based  their  derivations  on  the  work  of  Eckart.  The 


experimental  geometry  is  shown  in  Figure  4,  and  the  symbols  are  defined  as  follows: 

S  ,  =  highly  directional  acoustic  source 
d 

A  =  receiver 

X,  Y,  Z  =  orthogonal  Cartesian  coordinates  specified  by  unit 
vectors  i ,  j ,  q 

=  distance  from  the  source  to  the  origin  of  the  coordinates, 
which  is  positioned  in  the  center  of  the  insonified  area 
on  the  surface  S  [=  f  (x  ,  y ,  t)] 

R^  =  distance  from  the  origin  to  the  receiver 

r  =  xi  +  yj  +  zq  =  position  vector  from  the  origin  ro  an  elemental  area  dS 
on  the  surface,  which  has  a  unit  normal  n,  (positive 
down) 

R_  and  R.  =  distances  from  the  -«urce  5  ,  to  dS  and  from  dS  to  the 
S  A  d 

receiver  A,  respectively 
0  =  angle  of  incidence  of  the  acoustic  waves 

angle  of  scatter 

v  .j  =  angle  of  lateral  scattering  out  of  the  plane  of  incidence. 
Thesource  is  radiating  an  acoustic  field  <t>(x,  y,  z,  t),  which  is  given  as 

1 - 

<J)(x,  y,  z,  t)  =  11^  l(x,  y)  exp  i  jkR  -«  t  j  <  (12) 

R  V  2u  V 


where 


R  =  distance  from  S 
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II  =  source  power 
p  =  density  of  the  water 
C  =  speed  of  sound 

I(x,  y)  =  source  beam  pattern  (illumination  function) 
k  =  acoustic  wave  number 
jg  =  source  radian  frequency 
t  time 
i  . 

The  reradiated  acoustic  field  at  the  receiver  A,  which  is  produced  by 
reflection  and  scattering  from  the  rough  sea  surface  f(x,  y,  t)  is  given  by  the 
Helmholz  theorem^ 


<t>(A,  t)  =  --L  j  JdS 


<t> 


dn 


dn 


(13) 


where 


<t> 


ondi® 

dn 


=  values  of  <J>  and  its  normal  derivative  on  the  surface 


,  the  Green's 


The  value  of  <£ 


s  ond 


+  =  e«pi(k2.(RAC-o}/jRA0 

function  for  the  scattered  field 
kg  =  scattered  propagation  vector. 
d£> 


d  n 


can  be  determined  by  using  the  pressure  release 


condition  and  the  Kirchoff  approximation,  i.e., 


ad)  i 

1  d<&; 

=-d>.  and  — 

S  dn 

s  dn 

9 

S 

(U) 

where 


=  incident  acoustic  field  on  S 


=  normal  to  the  plane  tangent  to  the  point  of  evaluation  on  the  rough 
surface. 


Therefore, 


r)=-L  f  f dS-1  ($».*) 
4ir  Jr  J  3n 


Substituting  for  d>.  and  ^  and  partially  integrating,  as  in  Tolstoy  and  Clay, 
results  in 


<t>(A,  t) 


ik  E  f  (0)  exp  i  jk  (RAQ  +  R^ 


2*  RA0  RS0 


exp  (i  uQ  t) 


+co 

j  J dxdy  l(x,  y)  exp  i <|ax  +  by  +  cf(x,  y,  t)J  * 


where 


k=  I  k 


e-  m 
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f(9)  =  1 +  cos  £  cos  ^  ~  s'n  6  s'n  ^  cos 
cos  0  +  cos  ^ 

a  =  k(sin  0  -  sin  ^  cos  ) 
b  =  -k(sin  ^  sin  ) 
c  ~  -k(cos  0  +  cos  4>)  . 

in  the  derivation  of  Eq.  (16),  the  slopes  on  the  rough  surface  are  assumed  to  be 

small.  The  plane  wave  approximation  has  also  been  made,  since  differences  in 

direction  between  R_  and  R._  and  between  R.  and  R._  are  neglected. 

S  SO  A  AO 

Equation  (16)  cannot  be  applied  to  omnidirectional  or  wide-beam  sources, 
because  a,  b,  and  c  are  treated  as  constants  and,  hence,  the  insonified 
(illuminated)  area  on  the  surface  cannot  be  too  large.  In  addition,  the  curvature 
of  the  incident,  spherical  wave  fronts  must  be  very  small  so  that  the  plane  wave 
approximation  is  reasonably  good. 

g 

Clay,  however,  has  shown  that  this  theory  con  be  applied  to  the  case  of 
wide- beam  sources  by  subdividing  the  insonified  area  on  the  sea  surface  and 
applying  Eq.  (16)  to  each  subarea.  The  total  field  is  then  obtained  by  summing  the 
contributions  rrom  all  subareas  m  and  n.  The  scattering  geometry  for  the  long-crested 
sea  surface  is  shown  in  Figure  5,  where 

S  =  omnidirectional  source 

(0,  0)  =  originof  coordinates  halfway  between  the  source  and  the  receiver 
(**X,  D^),  (X,  D^)  =  source  and  receiver  coordinates,  'espectively. 

L  =  length  of  subsections 

(X  ,  0),  (X  ,  0)  =  centers  of  subsections  m  and  n,  respectively 
m  n 

Rc  ,  R_  =  distances  from  the  source  to  (X  ,  0)  and  (X  ,0),  respectively 
om  on  m  n 


R$m/  =  distances  from  the  source  to  (X  ,  0)  and  (Xn,  0),  respectively 

Ra  ,  Ra  =  distances  from  (X  ,  0)  and  (X  ,0),  respectively  to  the  receiver  A 
Am  An  m  n 

0,0=  angles  of  incidence  to  the  subsection  m  and  n,  respectively 
m  n 

^  ^  =  anglesof  scatter  from  subsection  mandn,  respectively,  to  the  receiver  A. 

When  Eq.  (16)  is  applied  to  individual  subsections  m  or  n  to  obtain  the  reradiated 

field,  <S>  (A,t)  or<t>  (A,t),  the  position  dependent  quantities  in  Eq.  (16)  must  be 
m  n 

appropriately  subscripted.  For  example,  when  a  isapplied  to  subsection  m,  a  becomes 

a  =  k(sin0  -  sin  t  ) . 
m  mm 

The  first  step  in  the  summing  procedure  is  to  separate  the  reradiated  field  into 

a  specular  (coherent)  component  and  a  scattered  (partially  coherent)  component, 

since  the  specular  contributions  add  coherently,  whereas  the  scattered  contribution? 
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are  only  partially  coherent  and  must  be  summed  separately.  This  means  that  the 
total  acoustic  field  <J>(A,  t)  is  given  by  the  sum  over  all  M  subsections: 


9 


(17) 


Mr  *i 

4>(A ,  t)  =  £  <t>  (A,  t)+  <J>$(A,  t) 

m=0  «-  c  J 


where 


<J>c(A,  t)  =  coherent  component  of  the  field 

<J>S(A,  t)  =  scattered  component  of  the  field. 

The  autocovariance  at  the  receiver  A  due  to  reflection  and  scattering  from  two 
subsections  m  and  n  is  given  by 


M  M 


/ <t> (A,  t)  <t>*  (A,  t  +  r)\  =  £  £  <4>  (A,  t)4>*  (A,  t  +  r  )\ 

X  X  m=On=0X  cm  cn  X 


M  M 


+  £  £  <  ^  (*'  *>  <*>*  (A,  t  +  r)^>,  (18) 

iifO  nm 


where 

*  denotes  the  complex  conjugate 
^  denotes  the  ensemble  average. 

This  result  follows  from  Eq.  (17)  because  the  mean  of  the  products  of  the  coherent 
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and  scattered  components  is  assumed  to  be  zero.  The  received  signal  spectrum 

can  now  be  obtained  by  substituting  Eq.  (16)  into  Eq.  (18)  and  Fourier  transforming 

the  result;  the  y  dependence  is  suppressed,  since  the  rough  surface  is  long- 

crested.  The  contribution  to  the  temporal  autocovariance  function  at  the  receiver 

A,  C  (T)  ,  from  ail  subsections  m  and  n,  is  the  ensemble  average  of  the 
mn 

product  of  <3^m(A/  t)  with  the  complex  conjugate  of  <J>^(A,  t)  : 


C(r  )  s  <«  (A,  f)  ®*  (A,  I  +  r  )  > 
inn  >  m  n 


k2E2f  »)!  (6) 
bi  n 

4*  Re  Ra  Rc  Ra 
Sm  Am  Sn  An 


exp 


?  {  MR; 


Sm  "  RSn  +  RAm  “  RAn*  “V*  } 


+oo 


7  /*“'*' 


“00 


I  W  !>')  exp  i  (o  x  -  a  x  ) 
m  n  m  n 


►<^expi  |  cmf^/  t)  -  cnt(*',  t+  r)|^> 


(19) 


where  the  primed  quantities  refer  to  the  nth  subsection.  The  time  average  has  been 
moved  inside  the  integral,  because  the  time  variation  of  the  surface  is  slow  com¬ 
pared  with  the  transmitted  signal. 

Since,  in  the  subsequent  analysis,  the  only  function  of  the  term  exp(-iu^r)  is 
to  locate  the  center  frequency  (carrier  frequency)  of  the  received  signal  spectrum, 
Uq  will  be  set  equal  to  zero,  and  the  calculated  spectrum  will  be  symmetric  about 
zero  frequency. 

Now,  by  performing  the  ensemble  average  on  Eq.  (19)  and  by  assuming  that  the 
distribution  of  surface  heights  :s  Gaussian,  we  see  that  the  autocovariance  becomes 

c  (  ,  .  fn<9>  ”P  !  {  k(RSm-  »Sn  *  «*■  -  *An>} 

4*2R<;  R.  Rc  R. 

Sm  Am  5n  An 


mn 


*•00 

*/  f* dx' 
-00 

exp  j-a2  [ 


I  (x)  I  (x')  exp  ?(a  x  -  ax) 
m  n  m  n 


2  2 

C  +  c  -  2c  C 
m  n  m 


'>] 


(20) 
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Since  (tie  illumination  functions  I  (x)  and  I  (x')  must  express  the  fact  that  the 

m  n 

acoustic  intensity  for  each  subsection  has  a  nonzero  value  only  within  the  region 


-L  to  L,  they  must  be  of  the  form 

=  =  <2,) 


where 


P  =  parameter  included  to  control  the  shape  of  the  function 
tb 

X  =  center  of  the  m  subsection 
m 

th 

X  =  center  of  the  n  subsection , 
n 

When  the  individual  subsection  illumination  (unctions  are  added  together,  the 
result  must  be  a  good  approximation  to  the  actual  source  beam  pattern  on  the  sea 
surface. 

A  change  of  variables  is  now  made  in  Eq.  (20)  in  order  to  express  the  auto¬ 
covariance  in  terms  of  space  lag  u. 


Let 


x  =  x"  +  X  -  - 
m  2 


x'=x"  +X  +- 
n  2 


(22) 


With  this  change  of  variables  in  the  integral,  the  limits  on  x7/  and  u  are  -co  to 
+  oo  and  the  Jacobian  of  the  transformation  is  one.  The  result  is 


C 

mn 


(0  = 


kVf  (8)  f  (9) 
m  n 

4tr^  R-  Ra  Rc  Ra 
jm  Am  Sn  An 


expi{k(RSm-RSn  +  RAm-RA„)} 
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Dropping  the  primes  and  factoring  out  the  constant  phase  term  results  in 


C  (r)  = 


k2E2f  (0)  f  (8) 

m  n 


4*  Rc  R.  R<.  R. 
5m  Am  5n  An 


CXP  5  |k(RSm'RSn  +  RAm’RAn)  +  {amXm‘°nXn)} 
(  im  on  Am  An  mm  n  n  l 


2  "i 

+  c  -  2c  C  p(u,  r  ) 
n  m  n  j 


(24) 


where 


Now  an  integration  over  x  can  be  effected  b/  factoring  and  interchanging  the 
order  of  integration  as  follows: 


C  (f)  =  G  exp  {-<r  (c 
(tin  mn  I  m 


-/)/ 


du|H(£  /  u)  exp(-iau) 


.exp|-2cmcn<.2[!-p(u,  r)]|  . 


Equation  (31)  is  analogous  to  Clay  and  Medwin's  Eq.  (16),  if  a  Gaussian  illumi¬ 
nation  function  is  substituted  for  I.  Witbin  the  limits  of  the  Kirchoff  approximation, 
Eq.  (31)  is  the  solution  for  the  covariance  of  the  received  acoustic  field  reradiated 
by  subsections  m  and  n  on  the  rough  surface.  Before  the  total  field  can  be 
evaluated  by  means  of  Eq.  (18),  the  covariance  must  be  separated  Into  coherent 
and  scattered  components.  This  is  done  by  making  the  :,$mail  roughness"  or  "long 
wavelength"  approximation^  and  retaining  only  the  first  two  terms  in  the  resulting 
power  series  expansion.  The  first  term  is  the  coherent  component,  whereas  the 
second  term  is  the  first-order  scattered  component.  For  larger  roughness,  one 
could  proceed  by  retaining  higher  order  terms  in  the  expansion.  The  small  rough¬ 
ness  approximation  is  made  as  follows: 


exP  !~2cm  V2  ['  -Ofr'  T>]  j  = 5 -2cm  V*!  e*P  j2 


Cm  Cnff2  Pk'r)i 
m  n  \ 


«exp(-2cm  cn<r2)  [^1  +  2Cffi  c.**  p(u,  r)J  ,(32) 


for  2  c  c  a  p(u,  r  )  <  <  1  . 
m  n  r  ' 
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where 

k2E2f  (0)  f  (0) 

G  =  m  -n 

mn  .  2 


A*  IL  R.  R*  Ra 

bm  Am  bn  An 


expi  jk(Rc  -Rc  +R.  -R,  )+(a  X  -a  X  )!  .  (27) 

r  |  Sm  Sn  Am  An  mm  n  n  j 


The  integral  over  x  can  be  defined  as 
+a> 

H(0,u)  =  /  dx  exp  (i  0x)  I 


(28) 


where 


0  =  (a  -  a  )  . 
m  n 


(29) 


After  this  substitution,  Eq.  (26)  becomes 


■  ao 


C  (r)  =  G  /  do 
mn  mn 


/ 


-co 


H(/3,  u)  exp(-iau) 


•  exp 


H 


r  2  +  2  , 

)  c  +  c  -  2c  c 
m  n  m 


np(x,  o]J 


(30) 


If  the  exponential  involving  p(u,  j- )  is  factored,  this  expression  can  be  written  as 


37 


Substituting  Eq.  (32)  into Eq.  (31)and  performing  appropriate  factoring  resuits  in 

+co 

Cmn(r)*Gmn  CXp|  [Cm  +  Cn j  |  J  du  H(d  '  u) 


-co 


+  2c  c  G  exp 
m  n  mn  r 


+0D 


/ 


-00 


du  H(0 ,  u)  exp(-iffu)  a  p(u,  r  ) 


(33) 


The  coherent  component  can  be  integrated  after  an  approprote  illumination 
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function  is  chosen.  As  in  Clay  and  Medwin,  the  Gaussian  illumination  function 
is  selected  in  order  to  eliminate  diffraction  sidelobes  resulting  from  the  source 
beam  pattern  and  to  facilitate  u  smooth  transition  between  subsections.  The 
function  H(&  ,  u)  is  evaluated  in  Appendix  A  ana  is  expressed  by 


H(0f  u)  =  PL (V2) 1/2  exp  ]-^2p2  L2/8  -  u2/2P2  L2  [  , 
which,  when  substituted  into  the  first  term  of  Eq.  (33),  results  in 

(A,  t)  $>*  (A.  t  +  t)  >  =  COH  J  C  (r)  =  P2  L2*  G 
\  cm  cn  s  j  mn  j 


mn 


exp 


2 

*  2  2] 

~a 

C  t  c 

i 

3 

3 

*1  -  (P2l2/4)(o^  «*> 

n  j  m  n 


(34) 


(35) 
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where  COH  jCffin(T)j  is  the  coherent  component  of  the  acoustic  autocovariance 
at  a  point  receiver  due  to  subsections  m  and  n  on  the  rough  surface.  From 
Eq.  (18),  the  total  coherent  component  is 


M  M 


-  (P2  L2/4)(a2  +  a2) 
m  n 


(36) 


The  coherent  component  is  thus  seen  to  be  independent  of  r  and,  hence,  does 
not  contribute  to  the  signal  fading.  The  spectrum  of  the  coherent  component, 
COH«{A(u)J  ,  is  obtained  by  talcing  the  Fourier  transform  of  C(r): 


r  MM 

COH  M(u)i  =  2P2  L2*2  £  EG  exp  j  -  o2  f(c2  +  c2)] 
'  '  L  m=0  n=0  mn  1  Lm  n  J 


-(P2  L2/4)(a2  +  a 
m 


^|]sw  , 


(37) 


which  shows  that  the  specular  (coherent)  component  is  an  impulse  at  the  trans¬ 
mitted  frequency,  for  conditions  of  iow  surface  roughness.  The  magnitude  of  this 
component  depends  upon  the  source  frequency,  the  experimental  geometry,  and 
the  mean-square  surface  roughness.  As  rhe  roughness  increases,  the  magnitude  of 
the  coherent  component  decreases. 

Application  of  Eq.  (37)  to  the  scattering  experiment  depicted  by  Figure  5 
shows  that  further  simplification  of  the  expression  for  the  coherent  component  is 
possible.  The  position-dependent  phase  term, 
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e*p  i  jk  [(RSm  +  RAm)  -  <RSn  +  RAn)]+  (am  Xm  -  an  X^)  j  , 

is  contained  in  G  of  Eq.  (37).  Since  there  is  typically  at  least  some  jitter  in 
mn 

the  hydrophone  position  in  measurements  at  sea,  the  R's  are  not  constant  but  vary 
randomly.  If  the  hydrophone  is  fixed  to  the  bottom,  the  R's  again  vary  randomly 
because  of  random  variations  in  the  ocean  stratification  that  change  effective  path 
lengths.  Hence,  it  is  desirable  to  perform  an  ensemble  average  on  this  term  over 
position.  This  average  is  not  performed  for  the  R's  in  the  denominator  because  the 
variation  with  change  in  R  is  very  slow  compared  with  the  exponent.  Now,  if 
the  position  dependent  variables  have  a  standard  deviation  such  that  the  phase 
terms  range  over  many  radians,  then 


^exp  i  kj(R 


'Sm  +  PAn)“(RSn  +  RAn)]' 


(a  X 
m  m 


0  if  m  ^  n 


1  if  m  =  n 


=  8  =  Kronecker  Deita, 

m,  n 


(38) 


because  the  phase  is  distributed  uniformly,  and  over  many  radians,  ihe  cross-product 
terms  nearly  cancel. 

Averaging  over  position-dependent  phase  and  substituting  from  Eq.  (27)  for 

G  yields 
mn  7 


COH 


(P2L2k2E2)/2 
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•z  z  fm<e>f„<e)^h2k+c„2)J 

m=0  n=0  '  L 


(P2L2/4)(o^o2)|ASmRAmRSnl<An|6W 


Grouping  the  m  ond  n  subscripted  expressions  results  in 

COh{a<,.»)}=  |{P2L2k2EV2 


•E  E 

m^O  n=0 


l  22  /D2 

*Xp\  Cm  ' 

L2  a2  /4)  f  (0) 
m  m 

RSm  RAm 

exp{-.2C2  -(P2  L2  a^/4) fn(0) 


RS.n  RAn 


6(u)  .  (40) 


A  A 

If  the  terms  in  braces  are  called  P  and  P  ,  respectively,  the  summation  term 

m  n 

in  Eq.  (39)  can  be  written  as  a  single  summation  over  i: 

MM  M  2 

E  E  PmP  =E  p.  . 

n&nZ  m  n  j~0  1 


Now,  from  the  geometry  of  Figure  5, 


fj(8)  = 


;) 
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and,  consequently, 


COHjA(u)|  =-P2L2k2E2 


1 


RSjRAi 


RAi  +  DlD2 


D,  Ra.  +  D 
1  A| 


•  exp 


-a2c?+  (P2L2a2)/4 


2 

S(u)  . 


(42) 


The  first-order  scattered  component  of  the  received  acoustic  signal  spectrum 
will  now  be  derived.  From  Eq.  (33),  the  scattered  component  of  the  autocovariance 
due  to  subsections  rn  and  n  is 


<*Sra  <*,.)*!„  (A,  .*')>  =SCAT  jcjr) 


=  2c  c  G  exp 
m  n  mn  r 


-  a 


k~m| 


+ao 

/ 


du  HO  ,  u)  exp(-iau)  <*  p(u,  t  )  . 


(43) 


-oo 


Now,  define  the  integral  as 


+co 

J  du  HO  ,  u)  exp(-iau)  <*2  p(u,  r  )  =  ff2  ©  (a,  0,  r  )  ,  (44) 

-00 

in  analogy  to  Parkins'  ^  procedure  for  the  case  of  single-area  scattering.  Then, 
substituting  Eq.  (44)  into  Eq.  (43),  we  obtain 
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SCAT  {C-JT)|=2cmc^  G_  exp  j-'2[c*  +  «*]  j  ff2«(a,  0,r)  .  (45) 


mn 


m  n  mn 


The  spectral  contribution  from  the  m^1  and  n*^  subsections,  SCAT  |Amn(o)j 
Is  determined  by  taking  the  Fourier  transform  of  Eq.  (45)  with  the  result  that 

SCAT  |  A  (u) | -  2c  c  G  exp  l-*2  [c2  +  c2]  j  ff2  W(«,  0,u)  ,  (46) 

1  mn  i  m  n  mn  {  L  m  nj  i 


where 


+oo 


o2  W( 


a  ,  3,  «)  2  J  drexp(iur)  e(a,  0,  r  )  . 


(47) 


-oo 
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The  long-crested,  linear,  Eulerian  model  of  random  gravity  waves  (Pierson,  '  and 

41 

Neumann  and  Pierson  )  is  employed  to  determine  p(u,  r)  in  order  to  evaluate 

2 

°  e(a,  0,  t)  . 

In  this  model,  lonf- crested  Airy  (small  amplitude)  waves  are  assumed  to  be 
propogating  parallel  to  an  acoustic:  transmitter- receiver  line.  The  sea  surface 
elevation  f(xj,  t-)  above  some  mean  level  f=0  is 


r( 


aj 


xj  -  nt  +  «||j2S(n)dn  , 


(48) 


where 


x  j  =  distance  along  the  source  receiver  line 


t.  =  time 


ft  =  surface  wave  radian  frequency 
g  =  acceleration  due  to  gravity 
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S(ft)  =  ocean  surface  frequency  spectrum 
«.  =  random  phase  . 

The  general  properties  of  this  model,  as  well  as  proofs  of  limit  theorems  per¬ 
taining  to  the  term  ^ 2S(S2)  dCl  ,  are  given  in  Pierson,^  Neumann  and  Piersonf  ^ 
42 

and  Kinsman.  The  n.^Jel  includes  the  following  features: 

1.  It  is  characterized  by  an  invariant  probabilistic  structure,  i.e.,  the 
statistical  properties  of  specific  samples  never  change. 

2.  The  waves  are  all  traveling  in  the  same  direction  (long- crested). 

3.  The  randomness  enters  the  model  through  the  phase  terms  t.  .  The  phase 

th 

of  each  (i  )  wave  is  between  0  and  2«  ,  with  ail  values  being  equally  likely, 
i.a.',  the  probability  of  occurrence  of  a  phase  angle  between  <t>  and  0  +  d0  is 


p(<J>  <t_  <0  +  d0)=—  .  (49) 

^  ‘  2w 

Thus,  each  surface  generated  by  Eq.  (48)  is  identical  except  for  the  phase. 

4.  The  model  is  strictly  valid  only  for  infinitely  deep  ocean  gravity  waves, 
where  the  frequency  n  and  the  wave  number  K  are  related  by  the  dispersion 
relation: 

ft2  =  9<  .  (50) 


From  Eq.  (48),  the  covariance  function  a  f(u,  r) 
the  ensemble  average  of  the  product  of  f(x.^ ,  t  j)  and 


is  then  formed  by  taking 

r(X|2>  «  follows: 


00 


f(xj2,fj2^  "  a  P(u'  r)  ~  J  6n  cos^-v  -  nT^$(fl)  f 


(51) 


where 


u=  x.}  -x.2 


r  =  f-i  -  t._ 

I1  |2 


fl=  radian  frequency 

S(n)  =  surface  wave -height  frequency  spectrum 


co 


■f 


*-=  i  6ns  (n) 

o 


p(u,  r)  *  autocorrelation  function. 

Equation  (51)  represents  the  model  as  it  will  be  used  in  further  applications  in 

this  study.  This  equation  represents  a  traveling  wave- 1  ike  correlation  function  in 

space- time.  If  bears  many  similarities  to  the  theoretical  correlation  functions 

20 

employed  in  the  acoustic  studies  of  Clay  and  Medwin  and  the  water  wave  studies 
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of  Clay  et  al  • 

Now,  Eq.  (51)  ;an  be  rewritten  as  the  sum  of  exponential  terms: 


(52) 
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When  this  result  is  substituted  into  Eq.  (44)  and  after  collecting  terms  on  u  ,  the 
result  is 


oo  +oo 


*2e( 


/dtfS(n)/doH«’'  U>“P  H-D-l  exp  j-i  SI  t 


-oo 


oo  +oo 


—  J  dfiS(li)  J  duH(0,  u)exp  -j^a+?.^u 


-00 


Now,  by  defining 


exp 


+oo 


ifl 


H  (&,  a  -  — ^  =  -1.  J  du  exp  -i  ^  u  H(0,  u) 


-oo 


+oo 


H^,a  +  3_^=_L  j  du  exp  -i^r  +  JL^u  H(0  ,  u) 


-00 


and  substituting  Eq.  (54)  into  Eq.  (53),  we  obtain 


oo 

a  e{a ,  ,  r)  =  ir  J d  Q  S(ft)  H  ^3,  «  -  — ^  exp(-i  fir) 

o  ^ 


00 


it  ^ dfiS(fi)  H  ,  a  ^  exp(ifli 


(53) 


(54) 


(55) 
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The  spectrum  function  a  W(a,  ,  u)  is  derived  by  taking  the  Fourier  transform  of 
Eq.  (55): 


+00  CO 

r^W{«/  0,  u)  =  ir  J  j  dr  dfi  exp  ji  (u -Q)r  |  S(R)  ^ 

“00  o 

+00  00 

+  ir  j'  ^  dr  dfl  exp  ji(o  +Q)r  $(R)  ,  a  +  —  ^\.  (56) 

-CD  O 


After  interchanging  the  order  of  integration  and  integrating  on  r,  the  result  is 


oo 


r2^*,  H,  «)  =  *  j  dflS(R)*2n6[u  -q]H 
o 


00 


+  *  j  d«s(n).2»6[u  +  n]  h 


/  2> 

[*,«-  — 

V  9> 

(s'“+  if)  ■ 


(57) 


Since  the  delta  function  "samples11  the  integral  at  fi  =  u  and  R  =  -u  , 
respectively,  it  follows  that 


r*W(«#  3,  u)  =  2»2  ,  a  -  — ^ 


S(u)  U  (m) 


+  2t2  H  ($,  a  +  —  ]  S(-u)  U  (-«)  , 
9  / 


(58) 


where 


U$j)  —  1  if  u  >  0,and  0  otherwise 
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U(-u)  =1  if  u  <  0/  and  0  otherwise  . 

Equation  (58)  can  be  evaluated  by  determining  the  functions  p,  .  This  can  be 
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done  by  a  method  due  to  Nuttall,  the  details  of  which  are  given  in  Appendix  A. 
The  result  is  that 


H 


4wP2L2  A 


-(PL0/2 


•  «St 


-(P  L  0  )/2  +  PL 


and 


2  2 

4*pV  A 


-(P  L  0)/2 


•'f') 


-(PL/8  )/2  + 


4  *; ) 


(59) 


where  <9  =  Fourier  transform  of  the  illumination  function.  For  the  Gaussian 
illumination  function,  Eq.  (59)  becomes  (Appendix  A) 


/  2\ 

r,  o  \2  1 

(p,  \  =P2L2exp 

^(p2L2/2) 

(«--  )  +*2/*\ 

V  9  / 

\  / 

'  9  '  J 

and 


(ff,«+S.\p,L2expj 

-(pV/^ 

(a  +  -2)  +d2/4 

'  9  /  1 

'  9  ' 

(60) 
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Now  after  substituting  Eq.  (60)  into  Eq.  (58),  applying  suitable  factoring,  and  re- 


After  substituting  Eq.  (61)  into  Eq.  (46)  and  summing  the  spectra!  contribu¬ 
tions  from  all  subsections,  the  scattered  received  signal  spectrum  SCAT  A(u)  is 
seen  to  be 


2  2  2  MM 

+  4,  P  l's<-«)U(-u)£  Y. 

m=0  n=0 


c  c  G  exp 
m  n  mn  r 


where  appropriate  factoring  has  been  performed. 
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Equation  (62),  Hie  scattered  acoustic  spectrum,  when  combined  with  Eq.  (42), 
the  coherent  component,  gives  the  complete  signal  spectrum.  The  scattered 
spectrum  is  symmetric  about  the  coherent  component,  with  an  up-Doppier  compo¬ 
nent  given  by  the  first  term  of  Eq.  (62)  and  a  down- Doppler  component  given  by 
the  second  term  of  Eq.  (62).  The  spectrum  of  the  scattered  component  thus  depends 
on  the  following: 

1 .  The  experimental  geometry 

2.  The  incident  radiation  wavelength 

3.  The  rms  surface  roughness 

4.  The  surface  height  frequency  spectrum. 

Frequencies  in  the  received  signal  spectrum  will  be  emphasized  when 


2  /a  +  a  \ 
u  .  /  m  n\ 

?  -i  )  ' 


because  the  functions 


exp  (  - 


2  2 
PL 


2  /a  +  a  \ 
u  /  m  nl 

g  \  2  / 


have  maxima  at  these  values.  These  are  the  resonance  relations  connecting  the 

frequency  and  geometry  of  the  experiment  with  the  component  waves  of  the  sea 

surface.  The  sums  of  the  direction  cosines,  (a  +  a  )/2,  behave  like  water  wave 

m  n 

numbers  in  determining  the  possible  scattered  frequencies  in  Hie  received  acoustic 

2  /o^  +  a  \ 

spectrum,  through  the  "dispersion  relation"  —  =  +  ( - 1  .  However, 

g  \  2  /  2  /a  +  a  \ 

since  the  exponential  functions  do  hove  nonzero  values  for  —  ^  +1 - PI  f 

9  \  2  / 
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small  contributions  will  enter  the  sum  from  tfiese  values  of  a  . 

Since  the  scattered  spectrum  predicted  by  this  equation  is  symmetric  about  the 
coherent  component,  it  is  necessary  to  consider  only  one  of  the  terms.  The  positive 
term  only  will  be  treated  in  the  work  that  follows.  This  term  will  be  defined  as 


Now,  when  Eq.  (64)  is  applied  to  the  scattering  experiment  depicted  in  Figure  5 
and  when  a  procedure  similar  to  that  for  the  coherent  component  is  followed,  i.  e. , 
Eqs.  (37)  to  (42),  the  result  is 


Aj{u)  =  (k4  EV  L2)  exp  -  (p2  LVAj2) 


M 

i=o 


F.(u) 


(65) 


where 


WD!V(IX 


.  exp  j-  O2  C2  -!P2L2a2)/4+P2l.2u2a  /2g  f  ,  (66) 

(  1  •'  i  \ 

Equation  (65)  represents  the  positive,  first-order  sideband  of  an  acoustic  signal  re- 
radiated  by  a  rough,  moving,  long-crested  sea  surface.  It  is  limited  to  cases  of  low 


surface  roughness  in  such  a  manner  that  the  term  2c  c  a  p(u,r)  must  be  <  1 

m  n 

in  Eq.  (32).  This  condition  is  reasonable  for  low  sea  state  underwater  sound 
experiments  using  low-frequency  acoustic  waves  at  large  angles  of  incidence 
relative  to  the  sea  surface.  Equation  (65)  will  be  a  better  representation  of  the 
scattering  from  the  swell  rather  than  a  representation  of  scattering  from  the, 
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locally  generated  sea,  since  the  swell  is  more  nearly  long- crested  than  is  the  sea. 

Equation  (65)  represents  the  scattered  acoustic  spectrum  as  the  product  of  the 
ocean  surface  wave-height  spectrum  and  a  weighting  function,  i.e. , 


Ajfc)*S(o)  B(u), 


(67) 


where 


Bfc)=Bk,  a.,  c.,  cr2)  =  (k4E2P2L2)eXpj-^2l2U4/2g2)jE  |f. 


fc>) 


,(68) 


which  depends  on  the  acoustic  frequency,  the  experimental  geometry ,and  the  mean- 
square  wave  height. 


3.3  THE  OCEAN  WAVE-HEIGHT  SPECTRUM 

IN  TERMS  OF  THE  ACOUSTIC  SCATTERED  SPECTRUM 


Since  the  acoustic,  first-order,  forward -scattered  spectrum  is  the  product  of 
the  ocean  wave-height  spectrum  5(u)  and  a  weighting  factor  B(u),  Eq.  (67),  it  is 
possible,  at  least  in  principle,  to  invert  Eq.  (67)  and  solve  for  $(»)  when  A(a)  is 
measured  and  B(u)  is  known: 


(69) 


S(u) =  . 

If  we  think  of  the  acoustic  source  and  receiver  setup  as  an  instrument  with  which 
to  measure  properties  of  the  ocean  surface,  then  B(u)  is  the  frequency  response 
function  of  the  "instrument. "  It  gives  us  the  frequency  band  on  the  surface  that 
can  be  measured  by  this  means,  as  well  as  the  gain  or  weighting  over  that  band. 

In  the  experiments  to  be  discussed  in  the  next  chapter,  measurements  of  the 
acoustic  forward- scattered  spectrum  and  the  ocean  wave-height  spectrum  made  in 
the  BIFI(Biock  Island- Fishers  Island)  acoustic  range  in  Block  Island  Souikj  will  be 
used  to  evaluate  the  usefulness  of  Eqs.  (6^  and  (69)  as  analytical  tools  for  study¬ 
ing  the  sea  surface. 
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4.0  FIELD  EXPERIMENTS  AND  RESULTS 
38  34 

Tolstoy  and  Clay  and  Fortum  have  summarized  the  results  of  the  most 
significant  model  tank  experiments  on  the  scattering  of  underwater  acoustic  waves 
from  rough  surfaces.  The  general  validity  of  both  Hie  physical  optics  and  the  wave 
expansion  theories  has  been  proved  for  low  surface  roughness,  even  under  conditions 
where  some  of  the  theoretical  assumptions  were  violated.  Thus,  sufficient  know¬ 
ledge  has  been  acquired  to  make  experiments  at  sea  both  practical  and  desirable, 
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especially  in  view  of  recent  advances  in  signal  processing  technology.  '  How¬ 
ever,  Marsh  and  Mellen,  ^  Parkins,  ^  and  Fortum*^  have  noted  that  few  bistatic 
measurements  of  surface  reflected  acoustic  signals  in  the  ocean  have  been  reported. 
In  addition,  in  many  of  the  sea  surface  scattering  experiments  described  in  the 
literature,  measurements  of  the  sea  surface  were  limited  to  visual  observations. 
Although  measurements  of  ocean  wave-height  and  surface- reflected  sound  pressure 
i9V‘r*  have  been  made,^'  a  greaf  need  still  exists  for  additional 

scattering  measurements  at  sea  under  conditions  varied  enough  to  test  the  applic¬ 
ability  of  existing  scattering  theory  to  ocean  wave  research. 

In  this  chapter,  the  results  from  several  acoustic  surface  scattering  experiments 
conducted  in  Block  Island  Sound  are  presented.  These  experiments  involved  sinu¬ 
soidal,  omnidirectional  acoustic  signals  transmitted  from  projectors  near  Block 
Island  and  received  by  means  of  two  widely  separated  hydrophones.  Ocean  wave 
height  was  monitored  simultaneously  by  a  resistance  wave  staff,  and  acoustic  and 
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ocean  wave-height  data  were  analyzed  for  spectral  content.  Results  from  the 
spectral  analysis  of  seven  of  these  experiments  are  presented  v/ith  particular 
attention  given  to  determining  relationships  between  the  ocean  wave-height 
spectra  and  the  first-order  sidebands  in  the  acoustic  spectra.  These  results  will  be 
used  to  evaluate  Eqs.  (67)  and  (69)  and  to  further  examine  the  potential  of  under¬ 
water  acoustics  as  a  too!  in  ocean  wind  wave  studies. 

4.1  AREA  OF  STUDY 

The  experimental  part  of  this  investigation  was  carried  out  in  the  BIFI  (Block 
Island  —Fishers  Island)  acoustic  range  in  Block  Island  Sound.  Block  Island  Sound 
(Figure  6)  is  situated  at  about  latitude  41°  101  N  and  is  bounded  by  Block  Island, 
Fishers  Island,  Long  Island,  and  the  Rhode  Island  shore  line.  It  is  shaped  like  a 
parallelepiped  and  has  a  surface  area  of  about  400  sq.  mi.  The  mean  depth  of 
Block  Island  Sound  is  40  m;  the  greatest  depth  is  about  100  m,  near  Fishers  Island. 

The  bottom  slopes  downward  from  Fishers  Island  and  the  mainland  toward 
Block  Island.  Near  Fishers  Island,  the  bottom  is  uneven  and  rocky  and  has  many 
depressions.  Near  Block  Island,  the  bottom  is  quite  flat  and  is  composed  of  sand 
and  mud. 

The  hydrographic  features  of  Block  Island  Sound  are  seasonal  and  thus  result  in 
(1)  negative  sound  speed  gradient  conditions  in  summer  and  early  fall  and  (2)  iso¬ 
velocity  or  slightly  positive  sound  speed  gradient  conditions  in  late  autumn,  winter, 
and  early  spring.  These  conditions  limit  surface  duct  acoustic  propagation  to  the 
cold  season.  The  results  of  a  veiocimeter  piobe  survey  made  in  January  1968  are 
shown  in  Figure  7.  This  section  exhibits  the  slightly  positive  sound  speed  gradients 
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that  are  typical  in  January  and  February.  The  horizontal  distribution  of  properties 
is  also  characterized  by  small  gradients  in  winter,  as  exemplified  by  Figure  8, 
which  presents  isotherms  contoured  from  aerial  rodiometer  data  taken  in  January. 

A  more  detailed  discussion  of  the  oceanographic  featuresof  Block  Island  Sound  is  g:ven 

51  52 

by  Williams,  Williams,  Azarovitz,  and  Lamoureaux,  and  Nalwalk,  Rathbun, 
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Rcbinson,  and  Riley.  Oceanographic  and  meteorological  observations  during  the 
time  of  the  measurements  are  given  in  Appendix  B. 

4.2  EXPERIMENTAL  SETUP  AND  INSTRUMENTATION 

The  facilities  of  the  BIFI  Range  in  Block  IslandSound  (Figure  6)  include  three 

bottom- mounted  receiving  hydrophones  at  Fishers  Islcmd,  three  bottom- mounted 

protectors  at  Block  Island,  and  various  types  cf  oceanographic  buoys.  A  complete 
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description  of  this  range  is  given  by  Hasse  and  Schumacher.  The  acoustic 
projectors  are  operated  from  a  manned  field  station  at  Block  Island,  and  receiving 
hydrophones  are  connected  by  cable  to  an  unmanned  field  station  at  Fishers  Island. 
Both  the  Block  Island  and  Fishers  Island  stations  are  linked  by  telephone  data  trans¬ 
mission  lines  to  the  BIFI  data  acquisition  and  analysis  laboratory  at  the  NUSC 
laboratory  in  New  London. 

Block  diagrams  showing  the  arrangement  of  the  transmitting  and  receiving 

equipment  used  in  the  experiments  are  shown  in  Figures  9  and  10.  Precise  signal 

frequency  control  was  provided  by  coherent  decade  frequency  synthesizers  at  the 

transmitting  and  receiving  locations.  These  instruments  are  equipped  with  tunable 

quartz-crystal -controlled  oscillators  with  resolution  to  nine  significant  figures. 

9 

They  have  rms  frequency  deviations  on  the  order  of  3  parts  in  10.  A  description 
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face  isotherms  contoured 
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Figure  9.  Block  diagram  of  the  signal  transmitting  equipment  on  Block  Island 
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of  the  instrumentation  used  in  this  study  is  given  in  Appendix  C. 

4.  3  PREUMiNARY  EXPERIMENTS 

A  number  of  preliminary  experiments  were  performed  to  determine  the  feasi¬ 
bility  of  making  a  detailed  study  of  the  effects  of  the  rough  sea  surface  on  sinu¬ 
soidal  signals  transmitted  across  the  B1FI  range.  Although  previous  studies  indicated 
that  surface  effects  would  dominate  the  short-term  signal  fading,  the  multipath  and 
multiple  reflection  propagation  conditions  in  the  BIF!  range  could  mask  the  surface 
scattering  effect. 

In  Hie  most  significant  of  these  experiments,  1702-Hz  tones  were  transmitted 

across  Block  Island  Sound,  from  Block  Islana  *o  Fishers  Island,  while  ocean  surface 

activity  was  monitored  by  an  accelerometer  wave  buoy  tethered  to  a  moored  buoy 

equipped  with  telemetry  apparatus.  Seventeen  runs  of  20-min  duration  each  were 

obtained  between  12  April  through  30  Mcy  1968.  The  following  conclusions  were 
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drawn  from  these  tests  : 

1.  Under  calm  sea  conditions,  no  frequency  spreading  occurred  above 
the  system  noise  level . 

2.  Under  negative  sound  speed  gradient  conditions  (late  May),  no  significant 
spreading  occurred,  even  though  the  sea  surface  was  quite  rough. 

3.  During  surface  duct  propagation  conditions  (rarly  April),  when  frequent 
interactions  between  the  acoustic  waves  and  the  sea  surface  occurred,  the  acoustic 
signal  spectrum  obtained  from  a  time  history  several  minutes  long  exhibited  first- 
order  sidebands  that  were  spread  in  frequency  from  the  carrier  by  an  amount  close 
to  the  peak  frequency  in  the  ocean  wind  wave  spectrum.  The  level  of  the  sidebands 


was  about  13  dB  below  the  carrier. 


Although  these  conclusions  could  be  drawn  from  the  1968  experiments,  it  was 
felt  that  greater  experimental  accuracy  could  be  obtained  in  the  ocean  wave 
measurements  by  using  a  resistance- hype  wave  staff.  In  addition,  acoustic  meas¬ 
urements  closer  to  the  transmitter  appeared  highly  desirable,  both  to  compare 
experimental  data  with  a  single- reflection  theoretical  model  and  to  examine 
qualitatively  the  effects  of  multiple  surface  reflections.  Hence,  a  series  of 
experiments  was  planned  for  the  winter  of  1969-70. 

4.4  ACOUSTIC  AND  OCEAN  WIND  WAVE 

MEASUREMENTS  IN  BLOCK  ISLAND  SOUND 

Sea  surface  elevation  and  underwater  surface- reflected  sinusoidal  signals  were 
obtained  simultaneously  at  both  short  (4.5  km)  and  long  (31  km)  ranges  from  the 
acoustic  projectors  at  two  different  frequencies  (127  Hz  and  17D2  Hz).  The  wave¬ 
lengths  of  these  signals  were,  on  the  average,  1 1 .4  m  and  0.86  m,  respectively. 
Both  acoustic  projectors  are  nearly  omnidirectional.  The  short-range  measurements 
were  taken  from  the  University  of  Connecticut's  research  vessel,  the  R/V  UCONN, 
a  20 -m  converted  Army  T-Boat. 

After  hydrographic  measurements  were  taken  at  Positions  B  and  C  (Figure  6),  a 
12-m  aluminum  spar  buoy  supporting  a  4.6-m  resistance  wave  staff  and  an  omni¬ 
directional,  bcoad-band  hydrophone  was  launched  at  Position  C.  The  buoy  was 
tethered  to  the  UCONN  by  means  of  a  76-m  mooring  and  signal  cable  with 
several  floats  attached.  The  configuration  of  the  buoy  system  is  shown  in  Figure  11. 
Its  resonant  period  of  oscillation  is  over  20  sec,  so  that  it  is  effectively  decoupled 
from  the  surface  wave  motion  prevalent  in  Block  island  Sound.  After  the  buoy 


Figure  1 1.  Spar  buoy  supporting  resistance  wave  staff  and  hydrophone 


Figure  12.  Spar  buoy  on  position  for  measurements  off  Block  Island 
(Photograph  courtesy  Prof.  D.  Paskausky.) 


The  wave  staff  and  hydrophone  were  then  calibrated,  and  a  20-min  recording  of 
both  the  wave  height  and  the  transmitted  acoustic  signal  was  made.  The  desire  for 
as  high  a  resolution  in  the  spectrum  as  possible  versus  the  length  of  time  that  the 
sea  surface  could  reasonably  be  considered  stationary  determined  the  selection  of 
the  20-min  time  history.  Prior  to  recording,  the  acoustic  signal  was  amplified. 
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filtered,  and  bond  shifted  to  8  Hz  to  faciiitote  analog -to -digital  conversion  and 
computer  processing.  The  wave-staff  signal  was  recorded  directly  on  magnetic 
tape  as  a  varying  dc  level.  During  the  later  experiments  in  this  series,  the  acoustic 
signal  was  transmitted  over  a  radio  frequency  link  and  was  received,  demodulated, 
bandshifted,  and  recorded  by  electronics  identical  to  those  used  for  the  underwater 
signal.  This  procedure  verified  that  none  of  the  sideband  structure  in  the 
acoustic  signal  spectrum  was  caused  by  instrumentation.  Upon  completion  of  the 
recording,  a  second  calibration  was  performed,  the  second  transmitter  was  switched 
on,  and  the  procedure  was  repeated.  Other  details  of  the  recording,  including 
block  diagrams  of  the  electronics,  are  given  in  Appendix  C. 

The  acoustic  signals  from  Block  Island  were  received  by  the  Fishers  Island 
hydrophone  and  transmitted  to  the  BIFI  laboratory  in  New  London,  where  they  were 
recorded  on  magnetic  tape  by  a  procedure  similar  to  that  on  the  UCONN.  After 
the  acoustic  and  ocean  wave  measurements  were  completed,  the  UCONN  steamed 
to  Position  A  and  another  hydrographic  cast  was  made. 


4.5  OPERATIONAL  RESULTS 

Fifteen  runs  were  made  in  Block  Island  Sound  between  29  December  1969  and 
30  January  1970.  During  these  runs,  simultaneous  measurements  of  wave  height 
and  sound  pressure  level  were  obtained,  and  the  acoustic  signal  from  the  Fishers 
Island  hydrophone  was  recorded  at  the  BIFI  laboratory  and  also  aboard  the  R/V 
UCONN.  Most  of  the  experiments  were  of  20  min  duration,  although  some  were 
shorter  because  of  operational  difficulties,  such  as  ships  moving  into  the  area. 
Although  all  the  experiments  produced  interesting  results,  this  discussion  is  limited 
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to  seven  cases,  which  are  sufficient  to  represent  a  variety  of  fair  weather  sea 
conditions  in  coastal  waters.  A  summary  of  the  operational  data  for  these 
experiments  is  given  in  Table  1. 

It  must  be  pointed  out  that  larger  sea  states  were,  of  course,  encountered 
during  storm  conditions,  but  when  such  conditions  occurred,  it  was  not  possible  to 
conduct  the  experiment  because  of  the  difficulty  of  holding  the  boat  in  position 
Gnd  of  launching  the  spar  buoy  from  a  pitching  deck.  In  addition,  for  wind  speeds 
greater  than  35  knots,  the  signal-to-noise  ratio  for  the  acoustic  signals  is  normally 
so  low  that  acoustic  reception  at  Fishers  Island  is  not  possible.  This  is  probably 
the  result  of  high  ambient  noise  levels  due  to  the  rough  sea,  as  well  as  to  the 
formation  of  a  bubble  layer  near  the  surface  of  the  water,  which  greatly  inoreases 
surface  loss  on  each  bounce.  Thus,  ail  the  data  presented  in  the  following  sections 
were  acquired  during  relatively  low  sea-state  conditions. 

4.6  DATA  PROCESSING  AND  ANALYSIS 

The  recorded  wave-staff  and  bandshifted  acoustic  data  were  signal-condi¬ 
tioned,  low-pass  filtered,  and  digitized  simultaneously  at  a  rate  of  64  Hz.  These 
procedures  were  followed  in  order  to  prevent  aliasing  of  the  computed  signal 
spectra.  Details  on  the  techniques  used  in  this  study  are  presented  in  Appendix  D. 

Variance  spectra  were  computed  by  means  of  the  Fast  Fourier  Transform  (FFT) 
algorithm  of  Cooley  and  Tukey  (Cooley,  Lewis,  and  Welch /.  Prior  to  spectral 
analysis,  the  dc  component  and  the  i  inear  trend  were  removed  from  the  data.  The 
time  histories  were  then  sectioned  into  overlapping  pieces,  which  were"  windowed” 
in  the  time  domain  to  minimize  spectral  leakage,  following  the  procedures  of 
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Welch  and  Bingham,  Godfrey,  and  Tukey.  The  FFfs  were  computed  (Arnold, 
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Nuttall,  Ferrie,  and  Carter  )  for  each  piece,  and  then  the  individual  estimates 
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were  averaged  to  achieve  statistical  stability  in  the  sense  of  Blackmon  and  Tukey. 

In  addition,  following  Jenkins  and  Watts, ^  window  carpentry  was  applied 
by  trying  different  windowing  functions,  and  window  closing  was  applied  to 
the  time  histories  by  computing  average  spectra  for  small  transform  sizes  (large 
bandwidth)  and  progressively  doubling  the  transform  size  until  maximum  resolution 
(smallest  bandwidth)  was  achieved.  The  increased  resolution  resulted  in  decreased 
statistical  stability. 

By  analyzing  many  samples  in  this  way,  it  became  clear  that  a  transform  size 
of  2,048  data  points  (resolution  of  0. 0225  Hz)  provided  the  best  trade  off  betveen 
resolution  and  stability.  These  considerations  are  dealt  with  in  greater  detail  in 
Appendix  D.  All  spectra  presented  in  this  section  are  resolved  to  0.0225  Hz. 

The  frequency  resolution  of  the  electronic  system  was  checked  by  recording 
and  processing  pure  tones  in  the  same  manner  as  the  data  signals.  A  resolution 
capability  of  better  than  0. 01  Hz  was  achieved .  A  field  check  was  obtained  by 
transmitting  127-Hz  and  1702-Hz  signals  from  Block  Island  to  the  UCONN  by 
means  of  a  radio  link.  The  signals  were  detected,  bandshifted,  and  recorded  by 
electronics  identical  to  the  water  path  system.  Figure  13A  shows  the  spectrum  of 
the  1702  Hz  signal,  and  Figure  13B  shows  the  spectrum  of  the  127  Hz  signal. 

These  tests  show  that  the  observed  variance  in  the  sidebands  of  the  acoustic  signal 
spectra  are  a  consequence  of  fluctuations  in  the  ocean  medium  and  are  not  due  to 
the  recording  electronics. 
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:igo.3  13  A.  Spectrum  level  of  1702-Hz  signal  transmitted  from  Block  Island 

to  VV  UCONN  via  RF  link 


RUNS  127  Hz 
R/V  UCONN 
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Figure  13  B.  Spectrum  level  of  127-Hz  signal  transmitted  from  Block  Island 

to  R/V  UCONN  via  RF  link 
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In  Hie  presentation  of  the  spectra  measured  from  the  UCONN,  the  ocean 
wave-height  spectrum  appears  directly  above  the  acoustic  signal  spectrum. 

The  variance  calculated  from  each  wave- height  record  is  presented  in  Table  2. 


Table  2 

VARIANCES  IN  THE  OCEAN  WAVE-HEIGHT  SPECTRA, 
RUNS  I  -  VII 


Run  Number 

Variance  (M^) 

1 

0.0591 

II 

0.0139 

III 

0.0032 

IV 

0.0021 

V 

0.0606 

VI 

0.0413 

VII 

0.0132 

The  center  frequency  of  the  acoustic  spectra  for  both  the  127- Hz  and  the  1702-Hz 
signals  is  8  Hz.  All  acoustic  signals  were  bandshifted  to  this  low  frequency  in 
order  to  facilitate  digital  analysis  of  the  20  min  data  samples  required  for  fine 
resolution  and  good  statistical  stability  (Appendix  D).  In  some  cases,  the  peak 
frequency  in  the  acoustic  spectra  are  not  exactly  at  8  Hz.  This  could  be  the 
result  of  three  factors:  (1)  failure  to  properly  tune  all  frequency  synthesizers 
(Appendix  C),  (2)  nonsrationary  processes  taking  place  in  the  ocean,  e.g. , 
accelerating  currents,  (3)  drift  in  the  position  of  the  spo  buoy  supporting  the 
hydrophone.  The  error  in  absolute  frequency  location  was  corrected  by  using  the 
peak  value  in  the  spectrum  as  the  reference  point  for  determining  Doppler  shifts 
(i.e.,  sideband  locations)  in  the  received  signal. 
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Since  the  sea  surface  modulates  the  transmitted  signal,  the  plots  of  the  acoustic 
spectra  are  double-sided,  i.e.,  sum  and  difference  frequencies  appear  about  the 
carrier  frequency.  On  the  other  hand,  the  ocean  wave-height  spectra  are  most 
properly  presented  as  single-sided,  positive- frequency  spectra.  However,  the 
acoustic  and  wave-height  spectra  can  be  compared  since  zero  frequency  in  the 
wave-height  spectrum  corresponds  to  the  bandshifted  carrier  frequency  (f^)  in  the 
acoustic  spectrum.  The  acoustic  up-Doppler  sideband  is  readily  compared  with  the 
wave-height  spectrum  by  aligning  frequency  f^  in  the  wave-height  spectrum  with 
frequency  in  the  acoustic  spectrum.  The  down- Doppler  component  is 

compared  by  taking  the  mirror  image  of  the  wave -height  spectrum  about  zero 
frequency  and  aligning  frequency  -f^  in  the  wave-height  spectrum  with  frequency 
fb  -  f  j  in  the  acoustic  spectrum.  In  the  following  section,  peak  locations  and 
slopes  in  the  acoustic  spectra  are  expressed  in  terms  of  the  magnitude  of  the 
frequency  displacement  away  from  the  carrier,  in  cases  of  asymmetry,  the  down- 
Doppler  component  is  used. 

4.  7  RESULTS  OF  THE  BLOCK  ISLAND  SOUND  EXPERIMENTS 

As  a  first  example,  consider  Run  I  (Table  1), which  was  obtained  near  the  end 
of  a  period  of  strong  northwesterly  winds  (Appendix  B).  The  wind-driven  sea 
surface  frequency  spectrum  was  single-peaked  with  maximum  energy  at  about 
0.25  Hz,  as  seen  in  Figure  14A.  No  swell  components  are  evident  in  this  spectrum. 
The  sea  surface  during  this  run  was  characterized  by  a  high  percentage  of  foam  and 
whitecaps.  The  acoustic  spectrum  shews  considerable  broadening  around  the  trans- 

19 

mitted  frequency  (the  location  of  the  specular  component).  Nuttall  and  Cron  refer 
to  this  as  "zeroeth  order  scattering."  A  distinct  down-Doppler  sideband  occurs 


Figure  14B.  Spectrum  level  of  acoustic  signal  at  Position  C  for  Run  I 
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33  dB  below  the  peak  level  of  the  carrier.  The  asymmetry  in  the  sidebands  is  prob¬ 
ably  due  to  multipath  interference. 

A  comparison  of  the  down-Doppler  scattered  spectrum  (defined  here  rather 
arbitrarily  as  between  7.  87  Hz  and  7.  25  Hz)  with  the  mirror  image  of  the  ocean 
wave-height  spectrum,  as  described  above,  shows  that  the  peaks  match  reasonably 
well.  The  acoustic  peak  is  0.016  Hz  (1  frequency  bin)  closer  to  the  reference  or 
"lero1*  frequency  than  is  the  wave-height  spectrum.  Since  the  two  spectra 
represent  estimates  from  random  processes,  this  small  difference  in  location  is 
probably  not  significant.  The  mean  slopes  of  the  spectra  are  nearly  equal  between 
0. 25  Hz  and  0. 5  Hz,  but  between  0. 5  Hz  and  0.  75  Hz  the  sea  surface  spectrum 
slopes  down  toward  higher  frequency  at  aboi't  24  dB/Hz,  whereas  the  acoustic 
spectrum  slopes  away  from  the  carrier  at  about  53  dB/Hz.  Frequencies  displaced 
by  more  than  +0.  75  Hz  from  the  carrier,  contain  no  information  since  the  noise 
floor  of  the  recording  and  processing  systems  has  been  reached. 

An  example  from  a  sea  surface  composed  of  a  swell  component  plus  residual 

wind  waves  is  given  by  Run  II.  In  this  case,  the  incoming  swell  was  running  almost 

counter  to  the  sea  (Table  1).  In  addition,  swell  waves  that  had  been  reflected  from 

the  beach  were  propagating  in  the  same  direction  as  the  wind  waves.  The  swell 

component  is  quite  possibly  the  result  of  a  storm  that  had  passed  off  the  Atlantic 
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coast  at  Norfolk  at  0700  hr  the  previous  day. 

The  spectra  at  Position  C  are  presented  in  Figure  15.  The  ocean  wave-height 
spectrum  exhibits  the  swell  peak  at  0.09  Hz  and  the  local  wind  wave  peak  at 
0.  34  Hz.  The  acoustic  spectrum  exhibits  a  greatly  broadened  central  peak,  i.e. , 
considerable  zero- order  scattering.  First-order  sidebands  occur  at  + .  08  Hz  from 
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Figure  15A.  Spectrum  level  of  ocean  wave  height  at  Position  ...  for  Run  II 
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Figure  15B.  Spectrum  level  of  acoustic  signal  at  Position  C  for  Run  II 
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the  carrier  of  about  6  dB  below  the  peak  level..  From  the  wave-height  spectrum, 
the  wind  wave  peak  is  5  dB  below  the  swell  peak,  whereas  the  corresponding  peaks 
in  the  acoustic  spectrum  differ  by  12  dB.  This  is  an  indication  that  the  low- 
frequency  wave  components  present  in  the  sea  surface  are  the  primary  forward 
scatterers  for  the  acoustic  waves.  Between  0.  5  Hz  and  0.  75  Hz,  the  wave- height 
spectrum  slopes  down  at  24  dB/Hz,  whereas  the  acoustic  spectrum  slopes  air 
18  dB/Hz.  Figure  16  shows  the  acoustic  spectrum  received  at  Fishers  Island  for 
Run  II.  This  spectrum  shows  a  narrower  specular  component  with  first-order 
sidebands  9  dB  and  11  dB  below  the  peak  carrier  level.  Between 0.  25 Hz  and  0.  5 Hz / 
the  slope  in  the  acoustic  spectrum  at  Position  C  is  about  18  dB/Hz,  whereas  the 
spectrum  at  Fishers  Island  slopes  at  36  dB/Hz.  Between  0.  5  Hz  and  0.  75  Hz,  both 
spectra  slope  at  18  dB/Hz. 
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Run  Hi  was  made  after  several  days  of  winds  from  the  west -northwest.  Before 
and  during  the  run,  the  sea  was  calming  down,  with  local  wind  waves  and  residual 
wind  waves  dominating  the  wave-height  spectrum  (Figure  17A),  which  peaks  at 
0.4  Hz.  Again  the  slope  between  0. 5  Hz  and  0.  75  Hz  (and  beyond  to  1  Hz)  is 
24  dB/Hz.  The  corresponding  acoustic  spectrum  (Figure  17B)  exhibits  broad  side¬ 
bands  having  a  peak  value  displaced  by  0. 4  Hz  from  the  carrier.  The  level  of  the 
sidebands  is  about  17dB  below  the  peak  carrier  level.  The  spectrum  at  Fishers 
Island  for  Run  III  (Figure  18)  exhibits  a  slightly  broadened  specular  component  but 
no  clearly  discernible  sidebands. 

Conditions  of  a  perfectly  flat  surface  would  have  been  desirable  as  a  reference 
for  these  experiments,  but  in  the  ocean  such  a  surface  never  exists.  The  closest 
approximation  to  a  completely  flat  surface  occurred  on  20  January  (Run  IV),  when 
a  very  slight  residual  swell  was  running  but  ripples  were  nearly  absent,  i.e. ,  the 
sea  was  becoming  glassy.  The  spectra  at  Position  C  are  shown  in  Figure  19.  The 
wave- height  spectrum  exhibits  the  swell  peak  at  about  0.  12  Hz.  Between  0. 5  Hz 
and  0.  75  Hz,  the  spectrum  slopes  down  at  32  dB/Hz,  which  is  considerably  steeper 
than  the  slopes  for  the  previous  three  runs.  This  is  the  result  of  the  near-glassy 
surface  conditions.  The  acoustic  spectrum  exhibits  very  weak  first-order  sidebands 
at  about  +0.  1  K  from  the  carrier.  The  level  of  the  sidebands  is  25  dB  below  the 
carrier.  Between  0.  5  and  0.  75  Hz,  the  spectrum  slopes  down  at  32  dB/Hz. 

The  spectrum  at  Fishers  Islcmd  (Figure  20)  exhibits  well-defined  first-order 
sidebands  at  about  +0.  12  Hz  from  the  carrier  and  16  dB  below  the  carrier  level. 

The  discrepancy  in  the  location  of  the  sideband  peak  in  the  acoustic  spectrum 
at  Position  C, as  compared  with  the  other  two  spectra,  could  be  due  to  the  fact  that 
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Figure  17A.  Spectrum  level  of  ocean  wave  height  at  Position  C  for  Run  III 
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Figure  17B.  Spectrum  level  of  acoustic  signal  at  Position  C  for  Run  III 
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Figure  18.  Spectrum  level  of  acoustic  signal  at  Fishers  Island  for  Run  III 


before  and  during  the  run  the  air  temperature  was  increasing  rapidly,  while  the 
water  surface  was  nearly  calm  (Appendix  B).  This  heating  could  result  in  a  warm 
surface  layer  that  would  refract  some  of  the  acoustic  energy  away  from  the  sea 
surface.  The  Fishers  lsla<xJ  spectrum  could  agree  with  the  location  of  the  swell 
peak  since  the  water  column  is  always  well  mixed  near  Fishers  Island  because  of 
high  current  speeds.  A  sample  of  the  actual  data  recorded  aboard  the  R/V  UCONN 
is  presented  in  Figure  21. 

Two  of  the  best  runs  in  this  series  were  made  on  26  January:  Run  V  at  1702  Hz 

and  Run  VI  at  127  Hz.  The  prevailing  weather  conditions  were  the  result  of  a 

stationary  front  that  had  developed  during  the  early  morning  hours  just  southeast  of 

Block  Island  and  that  covered  an  area  from  Atlantic  City,  N.  J. ,  to  a  low- 

63 

pressure  area  centered  in  the  Gulf  of  Maine.  The  spectra  at  Position  C  for 
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Figure  20.  Spectrum  level  of  ccoustic  signal  at  Fishers  Island  for  Run  IV 

Run  V  are  presented  i.i  Figure  22.  This  run  represents  the  heaviest  sea  conditions 
in  ihis  series  of  measurements.  Since  the  first-order  sidebands  at  +.12  Hz  from  the 
carrier  are  only  3  dB  below  the  carrier  level,  the  data  can  be  presented  cn  a 
linear  plot  (Figure  23).  There  is  good  agreement  between  the  location  of  the 
peak  in  the  wave-height  spectrum  and  the  first-order  sidebands  in  the  acoustic 
spectrum.  From  Figure  22,  the  slopes  of  both  acoustic  and  wave-height  spectra  are 
about  24  dB/Hz  between  0.  5and  0.75Hz.  In  Figure  23,  there  ore  smal  I  peaks  in  the 
acousticspectrumatabout0.25Hz;  these  peaks  are  AdBbelow  the  first-order  peaks. 
These  peaks  are  probably  clue  to  second -order  scattering  since  they  occur  at  twice 
the  frequency  of  the  first-order  peaks.  Evidence  of  a  small  amount  of  higher  order 

scatter  can  also  be  found  by  noticing  the  broadening  of  the  specular  component 
(although  this  broaden;^g  could  also  be  the  result  of  internal  waves).  The  spectrum 

at  Fishers  Island  is  shown  in  Figure  24.  The  specular  component  and  the  sidebands 
are  more  narrow  than  in  those  a*  Position  C,  but  the  possible  second-order  peaks 
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Figure  23A.  Spectra!  density  of  ocean  wave  height  at  Position  C  for  Run  V 
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Figure  23B.  Acoustic  intensity  per  unit  band  at  Position  C  for  Unit  V 
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Figure  24.  Spectrum  level  of  acoustic  signal  at  Fishers  Island  for  Run  V 

are  still  evident  at  6  dB  below  the  first-order  peaks.  Between  0.  5  Hz  and 0.75Hz, 
the  slope  of  this  spectrum  is  very  close  to  that  for  the  spectrum  at  Position  C,  i.  e. , 
24  dB/Hz.  Samples  of  the  data  recorded  during  Run  V  are  presented  in  Figure  25. 

Run  VI  was  made  approximately  1  hour  after  Run  V.  These  runs  illustrate 
the  difference  in  scattering  between  the  high-frequency  (1702  Hz)  and  the  low- 
frequency  (127  Hz)  acoustic  signals.  The  acoustic  spectrum  at  Position  C 
(Figure  26B)  exhibits  first-order  sidebands  at  + .  09  Hz  and  21  dB  below  the 
carrier,  as  compared  with  3  dB  for  the  1702-Hz  case.  The  sidebands  are  narrower 
and  exhibit  no  possible  second-order  peaks.  The  acoustic  sidebands  are  narrower 
than  the  wave-height  spectrum,  having  a  width  at  the  half-power  point  of  about 
0.06  Hz  as  compared  with  about  0. 1  Hz  for  the  wav-.-height  spectrum.  The  spec¬ 
trum  at  Fishers  Island  is  shown  in  Figure  27. 
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Figure  27.  Spectrum  level  of  acoustic  signal  at  Fishers  Island  for  Run  Vi 


Run  VII  took  place  at  Position  A  (Figure  6)  following  the  passage  of  a  cold 

64 

front  late  on  the  previous  day.  The  wind  direction  had  shifted  from  northeast  to 
northwest  prior  to  the  run  (Appendix  B).  Hie  wave-height  spectrum  (Figure  28) 
exhibits  considerable  energy  at  the  swell  band  (0.  08  -  0.  14  Hz),  plus  a  well- 
developed  wind  wave  spectrum  from  0.  24  Hz.  The  swell  direction  was  opposite  to 
that  of  the  local  sea  (Table  1).  The  very  low-frequency  peak  at  about  0.03  Hz  is 
due  not  to  the  wave  motion  but  rather  to  a  slow  oscillation  in  the  spar  buoy  induced 
by  the  extremely  strong  current  (greater  than  4  knots).  Low-frequency  sidebands 
about  23  dB  below  the  carrier  occur  in  the  acoustic  spectrum  at  close  to  +0.  10  Hz 
or  dose  to  the  swell  frequency.  The  wind  peak  at  0.27  Hz,  which  is  1  dB  higher 
than  the  swell  peak,  occurs  in  the  acoustic  spectrum  as  a  smill  peak  in  the  down- 
Doppler  sideband  8  dB  below  the  first-order  peak.  (This  small  peak  could  also 
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represent  second -order  scatter;  however,  second- order  scattering  at  this  frequency 
is  not  a  likely  possibility.)  The  Fishers  Island  spectrum  for  Run  VII,  which  is 
presented  in  Figure  29,  exhibits  most  of  the  same  features  as  the  spectrum  at 
Position  A. 


FREQUENCY  (Hr) 


Figure  29.  Spectrum  level  of  acoustic  signal  at  Fishers  island  for  Run  VII 


4.8  CORRECTION  ANALYSIS  OF  THE  SPECTRA 


A  quantitative  measure  of  the  degree  of  simile  ity  in  the  acoustic  scattered 
spectrum  A(f)  and  the  ocean  wave- height  spectrum  S(f)  can  be  obtained  by 
computing  the  coefficients  of  correlation  p  as  follows: 
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(70) 


where 

N  =  total  points  in  the  sum 
i  =  index  rf  summation 

k  =  starting  point  to  begin  the  computations  in  the  A  array.  (This  index 
is  used  to  align  the  scattered  part  of  the  acoustic  spectrum  with  the  ocean  wave- 
height  spectrum.) 

The  results  of  the  analysis  are  presented  in  Table  3. 


Table  3 

CORRELATION  OF  ACOUSTIC  AND  OCEAN  WAVE-HEIGHT 
SPECTRA  BETWEEN  FREQUENCY  INTERVALS 


Run  Number 

Frequency  Interval 

0-  0.5  Hz 

0.5-  1  Hz 

1 

0.97 

0.93 

II 

0.89 

0.96 

III 

0.97 

0.91 

IV 

0.75 

0.97 

V 

0.96 

0.94 

v\ 

0.93 

0.19 

VII 

0.63* 

0.73 

*  p  =  0.  98  between  0. 0  -  0.  15  Hz 
"p  =  0.  78  between  0.  15  -  0. 3  Hz 
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5.0  COMPARISON  OF  THEORY  AND  EXPERIMENT 

!n  this  chapter,  the  results  from  the  physical  optics,  long-crested  theoretical 
model,  as  expressed  by  Eqs.  (6 7)  and  (69),  are  compared  with  the  experimental 
measurements  presented  in  Chapter  4.  In  many  ways,  the  theory  is  incomplete,  in 
that  it  does  not  account  for  many  of  the  dyncmic  processes  taking  place  in  the 
ocean  or  for  the  multiple  reflections  from  the  boundaries  of  ihe  sea.  Nevertheless, 
the  acoustic  scattering  processes  predicted  by  the  theory  are  occurring  in  the  ocean, 
although  other  processes  not  taken  into  account  by  the  theory  are  also  taking  place. 

The  luxury  of  the  controlled  laboratory  experiment  is  not  usually  available  to 
the  geophysicist  in  the  field.  Uncontrolled  conditions  that  are  not  accounted  for  by 
theory  can  always  occur  during  the  experiment.  However,  the  lack  of  an  all- 
inclusive  theory  should  not  be  a  deterrent  to  the  application  of  a  simplified  theory  if 
useful  results  can  be  obtained.  The  results  presented  below  show  that  underwater 
acoustics  can  indeed  be  a  useful  remote  probe  of  the  ocean  surface,  if  the  proper 
precautions  are  taken  in  the  interpretation  of  the  measurements. 

5. 1  THE  SUITABILITY  OF  THE  PHYSICAL  OPTICS, 

LONG-CRESTED  MODEL  TO  THE 
BIFI  SCATTERING  GEOMETRY 

The  acoustic  experiments  in  Biock  Island  Sound  are  modeled  by  the  geometry 
shown  in  Figure  5.  This  idealized  model  neglects  acoustic  bottom  reflections,  as 
well  as  multiple  reflections  from  the  surface.  However,  under  surfcce-duct 
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propogaHon  conditions,  i.  e. ,  during  winter,  the  effects  of  the  bottom  ore  not  very 
significant  at  the  higher  acoustic  frequency  (1702  Hz).  ^  The  effects  of  multiple 
surface  reflections  of  the  acoustic  waves  propagating  from  Block  Island  to  Fishers 
island  were  reduced  by  taking  measurements  close  to  the  sound  source  at  Position  C 
(Figure  6).  Thus,  the  single-reflection  model  was  made  more  reasonable.  Exam¬ 
ine  m  of  the  data  in  Chapter  4,  however,  indicates  that  the  received  signal 
spectra  at  Position  C  and  at  Fishers  Island  are  quite  simiiarfjljhe  major  difference 


seems  to  be  that  the  bandwidth  of  the  received  signal  spectrumjjis  narrower  at 

Fishers  Island.  From  the  1 '  resonant  e' 1  relation  of  th  is  theory,*  /g  =  (a^  +  oJ/2  , 

higher  frequencies  in  the  received  signal  spectrum  correspond  to  steep  angles  of 

38 

scatter  that,  from  normal  mode  theory,  suffer  high  attenuation.  Thus,  a  suc¬ 
cession  of  surtace-boftom  reflections  would  attenuate  the  high-frequency 

31 

Doppler.  March  and  Kuo  have  shown  qualitatively  that  the  spectrum  of  a  signal 
scattered  many  times  from  the  sea  surface  should  net  differ  significantly  from  a 
signal  scattered  only  once. 

The  use  of  the  two-dimensional,  long-crested  scattering  mods!  cun  be  justi¬ 
fied  as  follows: 

1.  Fcr  a  source  and  receiver  at  long  range  and  shallow  depth,  the  Fresnel 

29 

zones  on  the  rurface  become  very  long  and  narrow  ;  hence,  the  most  important 
scattering  takes  place  along  a  relatively  narrow  strip  above  the  source-receiver 


2.  The  predominant  wave  directions  in  Block  Island  Sound  (west  to  west- 
northwest  for  wind  waves,  a.id  southeast  to  south  for  swell)  make  relatively  small 


angles  with  the  source- rece:ver  lines  (Figure  6).* 

3.  The  most  important  forward  scatterers  of  underwater  acoustic  waves  are  the 
long  swell  waves;  these  tend  to  be  more  nearly  long-crested  than  the  short  gravity 
waves. 

Perhaps  the  most  serious  criticism  that  can  be  leveled  against  the  application 

of  Eq.  (50)  to  the  Block  Island  Sound  data  is  that  the  long  swell  waves  in  the  area 

of  the  acoustic  range  are  net  '  *  infinitely  deep  water.  Strictly  speaking,  one  should 

42 

use  the  following  equation  in  place  of  Eq.  (50)  : 

=  gK  tanh  (Kh)  ,  (71) 


where 

h  =  water  depth  in  meters 
fi=2trf. 

However,  the  ensuing  mathematical  complications  would  greatly  reduce  the  useful¬ 
ness  of  the  model.  Therefore,  as  an  alternative  procedure,  Eq.  (50)  will  be  re¬ 
tained,  and  the  extent  of  error  examined.  Equction  (71)  was  used  to  obtain  the 
compute1-- generated  plot  shown  in.  Figure  30,  where  frequency  versus  wavelength 
curves  are  presented  for  water  depths  of  1C,  000  m  (in  effect,  infinitely  deep),  40m, 
30  mr  20  m,  10  m,  and  5  m. 

*Long  swells  from  the  aouthpjst  or  south  into  Block  Island  Sound  are  a'efocused 
by  the  Block  island  Channel  so  that  the  directions  of  propagation  make  reasonably 
small  angles  with  the  lines  "C"  -  "P"  and  "H"  -  ”Pn  in  Figure  6. 
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higure  30.  Frequency  versus  wovelength  for  ocean  £ 
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The  periods  of  swell  in  Block  Island  Sound  are  typically  on  the  order  of  10  sec 
(0.  1  Hz  in  frequency).  The  mean  water  depth  is  about  40  m.  Figure  30  shows  that 
a  10-sec  component  wave  (i.  e. ,  a  mathematical  Airy  wave,  not  a  real  wave)  has  o 
length  of  approximately  150  m  in  deep  water  and  a  length  of  144;  m  in  water  40  m 
deep.  More  importantly,  for  component  waves  of  a  given  length,  e.g. ,  150  m, 
there  is  a  2%  error  in  frequency  if  the  wave  is  assumed  to  be  in  deep  water  when  it 
is  really  in  water  40  m  deep.  This  error,  which  has  been  calculated  from  Figure  30, 
increases  as  water  depth  decreases  and  as  wavelength  increases.  For  wavelengths  of 
100  m  or  less,  there  is  a  very  small  error.  Figure  31  presents  a  computer- generated 
plot  of  the  phase  speed  squared  versus  the  water  wavelength.  There  is  a  3%  error 
in  phase  speed  if  a  150- m  wave  is  assumed  to  be  in  deep  water  when  it  really  is  in 
water  40  m  deep. 

Since  the  primary  purpose  of  this  study  is  to  investigate  a  genera!  method  rather 

than  to  obtain  the  best  possible  fit  to  a  particular  set  of  data,  the  deep  water  dis- 

2 

persion  relation  SI  =  gK  will  be  retained  (Figure  32)  with  the  knowledge  that 
there  is  some  error  present  in  the  modeling  of  the  low-frequency  swells,  the  extent 
of  which  can  be  estimated  from  Figures  30  and  31. 

2  I 

In  this  study,  the  small  roughness  approximation,  i.e. ,  exp  ^2cm  c^a  p(u,r)J 

2 

«  1  +  2c  c  a  p(u,  r  ),  was  utilized  because  all  the  data  cases  studied  were  with- 
m  n  r  ' 

in  its  limits  of  validity;  the  largest  error  encountered  (for  the  region  on  the  surface 
closest  to  the  source  and  receiver)  was  3%,  which  seems  reasonable  for  a  field 
experiment.  The  other  reason  for  using  the  small  roughness  approximation  is  that  it 
allows  the  ocean  surface  correlation  function  to  be  brought  out  of  the  exponent  and 


WAVELENGTH  (m) 

Figure  31 .  Water  wave  phase  speed  squared  versus  wavelength  for  ocean  gravity  (Airy)  waves  determined  from 

the  relation  v2=  (g/K)  tanh  (Kh) 
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into  a  linear  relationship  with  the  acoustic  covariance.  When  Fourier  transformed, 

this  relationship  permits  a  solution  for  the  ocean  wave- height  spectrum  directly  in 

terms  of  the  acoustic  scattered  spectrum,  which  is  the  primary  objective  of  this 

study.  An  alternative  approach  would  have  been  to  attempt  to  fit  an  assumed  form 

for  the  spectrum  (such  as  Pierson- Moskowitz  )  to  the  acoustic  data  and  minimize 

the  error  between  the  computed  and  measured  spectra  by  adjusting  the  parameters 

in  the  spectral  form  (i.e. ,  the  A  and  the  B  of  the  Pierson- Moskowitz^  spectrum). 

This  procedure  was  not  followed  here  because  a  direct  solution  was  desired  and 

because,  in  shallow  coastal  waters,  the  spectral  shape  is  not  readily  predicted  by 

any  of  the  standard  theoretical  ocean  wave  spectra  models. 

It  is  difficult  to  assess  the  importance  of  tidal  currents  and  internal  waves  in 

4 

Block  Island  Sound  on  hie  received  signal  spectrum.  As  Clark  and  Yarnall  have 
shown,  such  phenomena  produce  a  low-frequency  modulation  of  the  acoustic 
signal.  This  modulation  could  appear  as  a  broadening  of  the  coherent  component 
of  the  received  signal,  since  a  20- min  record  could  nut  resolve  these  low- 
frequency  components. 

5.2  THE  NUMERICAL  EVALUATION  OF  THE  ACOUSTIC 
AND  SEA  SURFACE  SPECTRA 

Equations  (67)  and  (69)  were  programmed  for  numerical  evaluation  on  the 

NUSC  UNIVAC  1108  digital  computer.  A  listing  of  the  program  is  given  in 

Appendix  E.  The  first  problem  to  be  overcome  in  the  evaluation  of  these  equations 

2 

is  the  determination  of  the  following  Darameters:  P  ,  the  illumination  function 
shape  control;  L,  the  length  of  the  subsections  on  the  surface;  and  M,  the 


number  of  subsections  used  in  the  evaluation. 


The  choice  for  P  is  fairly  straightforward  in  that  the  acoustic  "illumination" 
at  the  boundaries  between  subsections  should  not  be  very  different  from  that  at  the 
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centers  of  the  suosections.  If  the  value  of  P  is  not  chosen  properly,  a  "ripple" 

will  appear  in  the  illumination  and  thus  cause  unwanted  interference  effects.  When 

2 

P  =  1/2  ,  the  intensities  obtained  by  adding  the  overlapped  portions  of  the  segment 

illumination  functions  produce  a  reasonably  smooth  total  surface  illumination.  The 

29 

selection  of  L  and  M  is  more  difficult,  however.  Application  of  Beckmann's 

formulas  (9)  and  (10)  to  the  BIF1  geometry  yielded  a  Fresnel  zone  4,066  m  long. 

29 

Additional  discussions  in  Beckmann  and  S pi zzi chino  point  out  that,  for  weakly 

directional  antennas,  a  very  large  area  on  the  scattering  surface  is  involved  in 

8  30 

scattering  the  incident  radiation.  Clay  and  Medwin  and  Clay  show  that  in 
underwater  sound  experiments,  a  very  large  area  on  the  surface  is  involved  if  ail 
significant  scattered  contributions  are  considered.  On  the  basis  of  these  consider¬ 
ations,  the  scattering  area  for  this  model  was  made  as  large  as  possible  without 
violating  the  criterion  for  the  "small  roughness  approximation. It  was  arbitrar¬ 
ily  decided  to  keep  the  error  due  to  the  approximation  to  less  than  3%.  The  3% 
error  limitation  on  the  small  roughness  approximation  permits,  for  the  roughest 
surface  encountered  (Run  VI,  ^  ~  0.  06  m^)  ,  a  scattering  area  4,503  m  long. 


*An  easy  solution  to  the  determination  of  the  size  of  the  scattering  area  on  the 
surface  would  be  to  include  a!!  the  area  on  the  surface  between  the  source  and 
receiver.  However,  for  those  regions  on  the  surface  very  close  to  the  source  and 
receiver,  the  cm  in  the  small  roughness  approximation  are  large.  Hence,  for  the 
high- frequency  case,  this  criterion  would  be  violated. 
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For  the  case  of  the  1702-Hz  measurements,  this  region  comes  within  89  m  of  the 

source  and  receiver.  In  the  case  of  the  127-Hz  measurements,  it  is  possible  to 

consider  areas  much  closer  to  the  source  and  receiver  because  the  c  are  very 

m  ' 

small  as  a  result  of  the  small  wave  number  (k  =  0. 55).  However,  because  of  the 
long  wavelength  (1 1.4  m),  one  cannot  consider  scattering  areas  very  near  the 
source  and  receiver  without  violating  the  farfield  assumption.  Hence,  a  compro¬ 
mise  is  again  necessary  and  a  scattering  area  of  4,470  m  is  used  for  the  data  at  this 
wavelength. 

The  choice  of  L  is  limited  by  the  criterion  that  the  wavefront  be  essentially 

plane  over  the  subsection  length,  i.  e. ,  the  sums  of  the  direction  cosines  should  be 

20 

essentially  constant  over  the  length  L.  Clay  and  Medwin  used  the  additional 
assumption  that  the  size  of  the  subareas  should  be  large  compared  with  the  corre¬ 
lation  distance  on  the  water  surface.  This  assumption  was  necessary  in  order  to 
integrate  their  expression  for  the  scattered  covariance  using  a  theoretical  form  for 

the  surface  space-time  correlation  function.  Their  model  tank  experiments  did  not 

30 

quite  satisfy  this  assumption,  but  agreement  between  theory  and  experiment  was 
still  good.  In  the  present  study,  this  assumption  was  not  necessary  in  developing 
the  theory  because  the  scattering  equations,  Eqs.  (56)  and  (57),  could  be  integrated 
by  means  of  the  long-crested  surface  correlation  function,  Eq.  (51),  and  by  trans¬ 
forming  from  the  time  to  the  frequency  domain.  However,  by  averaging  Eqs.  (37) 
and  (62)  over  phase,  one  finds  that  the  cross-product  terms  between  subsections  are 
equal  to  zero.  This  is  tantamount  to  assuming  that  the  contributions  from  adjacent 
subsections  are  incoherent..  In  a  strict  sense,  the  reradiated  signals  from  adjacent 
subsections  can  be  incoherent  only  if  the  surface  correlation  distance  is  small  in 
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comparison  to  the  length  of  each  subsection  L.  However,  in  the  evaluation  of  the 
numerical  model,  Eq.  (67),  good  results  are  obtained  when  L  is  small  because,  in 
reality,  the  coherence  between  contributions  ftom  small  adjacent  subsections  is  low 
as  a  result  ■'f  “jitter"  in  the  receiver  location  and  of  nonstationary  conditions  in 
the  ocean. 

Marsh^has  argued  that,  in  numerical  models  of  this  type,  the  mesh  (subarea) 
size  should  be  made  very  small.  The  basis  for  this  agreement  is  that  the  double 
summation,  Eq.  (64),  in  three-dimensional  form  would  approach  an  integral  if  the 
mesh  size  is  made  arbitrarily  small.  If  the  integration  were  performed  over  all 
directions,  the  plane  waves  would  add  and  thus  synthesize  spherical  waves  radiating 
from  a  point  source  and  converging  on  a  point  receiver.  However,  in  this  particular 
numerical  model,  Eq.  (67),  there  is  a  lower  bound  on  the  size  of  L,-  which  is 
determined  by  the  requirement  that  Hie  phase  terms  in  Eq.  (64)  average  to  zero  to  a 
good  order  of  approximation.  This  means  that  L  must  not  be  so  small  that  contri¬ 
butions  from  adjacent  subsections  are  coherent.  In  order  to  prevent  this  from 
happening,  the  lower  bound  on  L  has  been  set  by  requiring  that  the  length  of  the 
subsections  be  equal  to  several  radiation  wavelengths. 

The  criterion  that  the  sums  of  the  direction  cosines  be  relatively  constant  over 

each  subsection  was  interpreted  to  mean  that  the  change  in  the  value  of  a  or  c 

m  m 

from  one  subsection  to  the  next  should  be  less  than  10%  .  This  constraint  limits  the 
size  of  L  to  about  5  m  at  1702  Hz.  For  the  1702-Hz  data,  a  value  of  L  =  3  m  was 
adopted,  which  is  almost  four  times  the  acoustic  v/avelength. 

For  the  127-Hz  case,  the  wavelength  is  11.4  m  and  the  criterion  that  the  a 

m 

differ  from  section  to  section  by  less  than  10%  cannot  be  met  without  making  L  the 
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size  of  a  wavelength  or  less,  which  introduces  some  conceptual  problems  in  apply¬ 
ing  the  theory.  In  an  attempt  to  keep  the  ratio  of  section  length  to  wavelength 
constant  for  both  data  cases,  a  value  of  L  =  30  m  was  selected  for  the  127-Hz 
data,  although  the  applicability  of  the  theory  to  this  case  is  marginal.  Values  of 
L  equal  to  15,  20,  25,  30,  and  50  m  were  tried  for  this  cu  e  without  much  change 
in  result.  For  the  1702- Hz  case,  several  values  of  L  were  also  tried  up  to  50  m. 
The  numerical  results  at  3702  Hz  were  reasonably  consistent  for  values  of  L  below 
10  m. 

The  numerical  results  for  the  1702-Hz  ar.d  127-Hz  data  will  now  be  presented. 

In  this  chapter,  all  formulas  expressed  in  Chapter  3  in  terms  of  radian  frequency 

have  been  converted  to  frequency  f.  The  mirror  image  of  the  down -Doppler 

measured  acoustic  spectrum  is  used  in  the  calculations. 

The  weighting  or  "transfer"  function  B(f),  in  Eqs.  (67)  and  (69),  depends  on 

2  2 

the  mean-square  wave  height  o  .  In  this  study,  a  was  measured  directly  by 
wave  staff,  although  in  principle  it  could  be  determined  from  the  acoustic  meas¬ 
urements  from  the  coherent  component  by  means  of  Eq.  (42).  Through  an  iteration 
2 

procedure,  <r  could  be  adjusted  until  the  error  between  the  computed  and 
measured  coherent  components  is  a  minimum  in  the  least-squares  sense. 

The  function  B(f)  normalized  to  the  peak  value  is  shown  in  Figure  33.  This 
particular  B  function  has  been  computed  from  the  data  for  Run  III,  but  no  detect¬ 
able  differences  in  the  plots  can  be  seen  in  the  B’s  computed  from  the  other  runs. 
The  B  function  has  also  been  estimated  by  dividing  the  measured  acoustic  spectrum 
A^(f)  by  the  measured  wove- height  spectrum  S(f)  .  Thus,  from  Eq.  (67),  B (f)  is 
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Figure  33.  Normalized  B  function  for  acoustic  wave  number  k  =7. 38/m; 

Block  Island  Sound.  This  functic  o  was  computed  using  data  from  Run  III 

B(0  =  Aj(f)/S(f)  .  (72) 

Since  the  individual  measured  B's  are  quotients  of  two  random  processes  and 
thus  tend  to  be  very  noisy,  the  average  B  has  been  obtained  for  each  wavelength. 
The  B  for  the  1702-Hz  case  is  presented  in  Figure  34,  where  the  function  begins 
at  about  0.  1  Hz  because,  below  0.  1  Hz,  energy  from  the  broadened  "specular" 
component  can  be  found  which  invalidates  this  calculation.  The  general  form  of 
the  measured  and  computed  B's  is  similar;  however,  above  0.  3  Hz  the  computed  B 
falls  off  much  more  rapidly  than  does  the  measured  B. 

The  first  example  to  be  considered  is  Run  II.  The  acoustic  scattered  spectra, 
normalized  to  the  peak  values,  are  shown  in  Figure  35A.  The  measured  spectrum 
is  shown  only  above  0.05  Hz  because  the  coherent  component  is  encountered  below 
that  frequency.  The  computed  spectrum  has  been  evaluated  from  Eq.  (67)., 


Figure  34.  Normalized  average  measured  B  function  for  acoustic  wave 
number  k  =  7. 3E/m;  Block  Island  Sound 

The  normalized  ocean  wave- height  spectra,  measured  and  computed  by  means 
of  Eq.  (69),  are  compared  in  Figure  35B.  The  computed  ocean  wave-height  spec¬ 
trum  is  bandlimited  because  the  B  function  approaches  zero  as  0. 5  Hz  is  approached 
(Figure  33).  Hence,  in  this  case  and  in  the  following  cases,  the  A^(tyB(f)  compu¬ 
tation  is  truncated  at  the  point  where  it  starts  to  increase  rapidly.  The  agreement 
between  peak  values  is  generally  good.  Since  both  the  measured  curves  and  the 
computed  curves  are  estimates  from  random  processes,  the  alignment  is  not  expected 
to  be  perfect.  The  measured  acoustic  spectrum  does  not  reflect  the  small  peak  at 
0. 34  Hz  in  the  computed  acoustic  spectrum.  Consequently,  the  peak  at  0. 34  Hz  in 
the  wave-height  spectrum  is  not  detected  in  the  computed  spectrum. 

The  spectra  for  Run  ill  are  shown  in  Figure  36.  There  is  generally  good  agree¬ 
ment  out  to  0.44  Hz.  This  is  an  interesting  case  in  that  the  ocean  wave-height 
spectrum  pecks  at  an  unusually  high  frequency  for  Block  Islnnd  Sound. 
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Figure  35A.  Normalized  acoustic  spectra  for  Run  II,  1702  Hz 


Figure  358.  Normalized  ocean  wave-height  spectra  for  Run  li 
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The  results  of  Run  SV  are  given  in  Figure  37.  Agreement  between  theory  and 
experiment  is  good  between  0.  1  Hz  and  0. 35  Hz.  However,  the  double  peak 
below  0.  1  Hz  In  the  measured  acoustic  spectrum  may  be  due  to  either  statistical 
nof<e  or  low-frequency  modulation  by  infragravity  waves.  This  lower  peak  is 
falsely  predicted  in  the  computed  ocean  wave-height  spectrum. 

A  relatively  rough  surface  case  is  shown  in  Figure  38.  (Before  this  run  began, 
a  number  of  improvements  had  been  made  in  the  electronics  to  increase  the  signal - 
to- noise  ratio  (Appendix  B)  of  the  measured  data.)  There  is  good  agreement 
between  the  ocean  wave-height  spectra  out  to  0. 46  Hz,  where  the  Aj(f)/B(f)  compu¬ 
tation  is  starting  to  increase  rapidly. 

The  B  function  for  the  127-Hzdato  ?s  shown  inFigure  39.  Itis  much  narrower  than 
the  1702-Hz  B  function  and  indicates  that  only  ocean  waves  below  0.2  Hz  are  effec¬ 
tive  in  modulating  the  acoustic  signal.  The  narrowness  of  the  B  function  at  127  Hz  im¬ 
plies  that  information  about  the  sea  surface  obtained  by  acoustic  scatter  measurements 
is  bandiimited  to  0.2  Hz  for  the  BIFI  geometry.  The  measured  average  B  function  is 
shown  in  Fig.  40.  Again,  the  computed  B  falls  off  more  rapidly  than  the  measured  B. 

The  first  127-Hz  data  case  to  be  considered  it  Run  1  (Figure  41).  Agreement 
between  theory  and  experiment  is  very  poor.  During  this  run,  the  sea  surface  was 
characterized  by  a  high  incidence  of  whitecaps,  foam,  flying  spray,  and  turbulence, 
which  may  have  acted  in  concert  to  invalidate  the  assumptions  of  the  theory. 

The  situation  for  Run  VI  is  somewhat  better  (Figure  42).  The  measured  acoustic 
spectrum  is  severely  bandiimited  as  compared  with  the  ocean  wave- height  spectrum. 
This  result,  which  would  be  expected  from  the  B  function  (Figure  39),  limits  the 
usefulness  of  this  frequency  in  measuring  the  ocean  wave-height  spectrum,  as  can 
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Figure  37A.  Normalized  acoustic  spectra  for  Run  IV,  1702  Hz 
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Figure  37B.  Normalized  ocean  wave-height  spectra  for  Run  IV 
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Figure  38A.  Normalized  acoustic  spectra  for  Run  V,  1702  Hz 
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Figure  388.  Normalized  ocean  wave-height  spectra  for  Run  V 
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Figure  39.  Normalized  B  function  for  acoustic  wave  number 
k  =  0.55/m;  Block  island  Sound.  This  function  was  computed 
using  data  from  Run  VI 

be  seen  in  Figure  41 B.  The  low-frequency  end  of  the  peak  is  well  reproduced,  but 
no  information  above  0.  18  Hz  can  be  recovered.  This  effect  is  more  dramatically 
illustrated  in  Run  VII  (Figure  43),  where  the  false  peak  in  the  ocean  wave-height 
spectrum  at  0.03  Hz  in  the  measured  spectrum  has  been  filtered  out.  The  low- 
frequency  swell  peak  at  0.  1  Hz  is  predicted  almost  perfectly  by  the  acoustic 
method,  but  the  wind  wave  peak  at  0.  25  Hz  is  completely  absent. 

5.3  DISCUSSION 

There  is  better  agreement  between  theory  and  experiment  at  1702  Hz  than  at 
127  Hz.  This  result  is  to  be  expected  since  the  applicability  of  the  physical  optics 
theory  to  the  127-Hz  BIF1  scattering  experiments  is  marginal.  At  127 Hz,  the 
acoustic  wavelength  of  greater  than  11  m  is  about  one- fourth  of  the  depth  of  the 
water,  and  bottom  effects  become  very  significant  even  under  positive  gradient 
conditions.  For  this  low  frequency,  a  normal  mode  model  in  which  the  scattering 
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Figure  40.  Normalized  average  measured  B  function  for  acoustic 
wave  number  k  =  0.55/m;  Block  Island  Sound 

effects  enter  the  problem  through  the  boundary  conditions  Ytould  most  probably  yield 
better  results. 

However,  from  both  the  measured  and  computed  weighting  functions  B(f) 
(Figures  39  and  40),  it  seems  clear  that  this  low  acoustic  frequency  discriminates 
strongly  against  the  high-frequency  components  on  the  sea  surface,  as  would  be 
expected  from  the  Rayleigh  criterion.  Another  factor  that  significantly  reduces 
the  usefulness  of  this  low-frequency  sound  in  wind  wave  studies  is  that  the  strength 
of  the  surface  wave  modulation  is  frequency  dependent  (Eq.  68)  and  decreases  with 
increasing  acoustic  wavelength.  This  effect  is  graphically  illustrated  in  Figures 
22B  and  26B.  For  a  very  similar  sea  surface  condition,  the  scattered  acoustic  spec¬ 
trum  was  3  dB  below  the  specular  component  at  1702  Hz  but  21  dB  below  the 
specular  component  at  127  Hz.  For  processing  systems  with  limited  dynamic  range, 
the  127-Hz  sidebands  would  be  in  the  noise  level  of  the  system  and,  hence,  of  no 
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Figure  41A.  Normalized  acoustic  spectra  for  Run  I,  127  Hz 
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Figure  41 B.  Normalized  ocean  wave-height  spectra  for  Run  I 
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Figure  42A.:  Normalized  acoustic  spectre  for  Run  VI,  127 Hz 
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Comparison  of  the  computed  and  measured  6  functions  for  both  1702-Hz  and 
127-Hz  signals  reveal  that  the  computed  B  function  falls  off  with  frequency  more 
rapidly  than  does  the  measured  B,  The  narrowness  of  the  computed  B  function  is 
most  probably  due  to  the  fact  that  a  long-crested  surface  wave  model  was  used  to 
predict  the  field  scattered  from  a  sea  surface  that  is,  in  reality,  short- crested. 

The  long-crested  wave  model  Eq.  (51)  has  interpreted  all  the  waves  "seen”  by  the 
wave  staff  to  be  propagating  in  a  direction  parallel  to  the  source -receiver  line. 
Hence,  wave  trains  propagating  at  some  angle  to  the  source- receiver  line  would 
be  modeled  as  having  greater  wavelengths  than  they  actually  have  (greater  by 
A/cos  0  ,  where  0  >s  the  angle  of  the  waves  with  the  source- receiver  line  and 
A  is  the  surface  wavelength).  Thus,  it  follows  that  the  surface  long-crested  co- 
variance  function  Eq.  (51)  overemphasizes  the  long  wavelengths  at  the  expense 
of  the  short  wavelengths.  Therefore,  the  long-crested  scattering  model  does  not 
adequately  represent  the  effects  of  the  higher  frequency  surface  waves.  Hence, 
since  the  B(f)  function  is  calculated  by  assuming  that  ai!  the  waves  are  propagating 
in  a  direction  parallel  to  the  source-receiver  line  when  they  really  are  not,  it  fol¬ 
lows  that  B(f)  will  be  deficient  at  high  frequencies  to  an  extent  controlled  by  the 
deviation  of  the  real  sea  surface  from  the  ideal  long-crested  sea  surface. 

In  addition,  the  long- crested  scattering  model  does  not  td<e  into  account 
scattered  contributions  from  outside  the  plane  of  incidence.  These  contributions  may 
be  characterized  by  large  ocartering  eagles,  which  correspond  to  high  frequencies. 


On  the  basis  of  these  considerations,  it  is  highly  probable  that,-  when  this 
physical  optics  mode!  is  generalized  to  include  the  short-crested  sea  surface  model 
(Pierson^),  the  slope  in  the  B  function  (Figures  33  and  39)  will  be  reduced.  The 
deficiency  of  the  model  at  high  frequencies  offers  a  possible  explanation  for  the 
poor  results  of  Run  I  and  the  reasonably  good  results  of  Run  VI.  This  effect  is 
illustrated  by  comparing  the  127- H2  acoustic  spectra  for  Runs  I  and  VI  (Figures  41A 
and  42A).  During  Run  1,  the  sea  surface  was  wind  driven  and  consisted  primarily 
of  short,  choppy  waves.  During  Run  VI,  the  saa  surface  was  dominated  by  incoming 
swells.  Since  a  much  better  approximation  of  a  long-cresred  sea  is  provided  by 
swells  than  by  short  wind  waves,  the  results  from  Run  VI  are  significantly  better  than 
those  for  Run  I.  Thus,  the  agreement  between  measured  and  predicted  acoustic 
spectra  is  partially  a  measure  of  how  nearly  long- crested  is  the  sea  surface. 

This  effect  is  not  so  pronounced  for  the  1 702—  Hz  experiments  since  the  shorter 
acoustic  wavelength  used  gives  greater  weight  to  the  short  gravity  waves  on  the  sea 
surface.  Hence,  good  results  were  obtained  for  Rur,  Hi  e^en  though  the  wave- 
height  spectral  peak  occurred  at  0.4  Hz. 

When  the  physical  optics  short-crested  model  is  developed,  the  B 
function  (Figures  33  end  39)  will  probably  approach  unity  for  the  case  of 
an  omnidirectional  source  and  receiver  located  far  from  land  boundaries. 
Therefore,  it  is  reasonoble  to  conclude  that,  when  all  angles  of  incidence 
and  reflection  are  taken  into  account  and  a  short-crested  model  for  the 
sea  surface  is  employed,  the  generalized  physical  optics  model  will  be  in 
agreement  with  Marsh's  prediction  that  the  first-order  sidebands  are  the 


uniformly  weighted  sea  surface  height  spectrum.  Hence,  the  use  of  an  omnidirec¬ 
tional  or  wide- beam  source,  combined  with  the  generalized  physical  optics  theory, 
greatly  increases  our  ability  to  obtain  sea  surface  information  by  acoustic  means. 

5.4  PHYSICAL  OPTICS  THEORY  AS  AN  AID 
TO  OCEAN  WAVE  STUDIES 

Parkins'  single-area  scatter  theory^  predicts  that  the  reradiated  spectrum 
consists  of  a  line  at  the  carrier  frequency  plus  two  delta  functions  that  depend  on 
the  acoustic  wavelength,  the  surface  wave-height  spectrum,  and  the  angle  of 
incidence  and  observation.  Hence,  for  the  low  roughness  case.  Parkins'  theory 
could  yield  information  about  only  one  frequency  component  in  the  sea  surface. 

The  generalized  physical  optics  theory,  as  presented  here  for  the  omnidirectional 
source,  considers  contributions  from  many  subsections  where  each  subsection  corre¬ 
sponds  to  a  Parkins'  type  model.  When  the  contributions  from  all  the  subsections 

20 

are  summed  by  the  method  r, f  Clay  and  Medwin,  the  result  is  not  three  delta 
functions  but  a  spectral  band  with  contributions  coming  from  each  elementary 
section.  Each  section  contains  information  about  a  particular  frequency  component 
of  the  sea  surface  so  that,  in  theory,  a  wide  band  of  surface  frequencies  can  be 
measured. 

The  nature  of  the  transfer  function  B  of  the  acoustic  experiment  is  that  of  a 
low-pass  filter  that  is  relatively  flat  out  to  a  certain  point  and  then  begins  to  fall 
off.  It  appears  that  the  rate  of  fall  predicted  by  the  physical  optics  long-crestec! 
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model  is  too  high  urd  would  be  decreased  by  employing  the  short-crested  model .  ' 

These  considerations  from  the  physical  optics  method  indicate  that  the  most 
practical  way  to  measure  the  wave-height  spectrum  acoustically  is  to  spectrally 
analyze  the  surface  reradiated  signal  from  an  omnidirectional  source.  If  the 
surface  is  not  too  rough,  the  distribution  of  variance  in  first-order  sidebands  in  the 
acoustic  spectrum  will  be  similar  to  the  sea  surface  spectrum  out  toagiven  point  on  the 
frequency  axis  when  the  weighting  function  begins  to  fall.  When  it  is  necr  zero, 
no  information  about  the  wave-height  spectrum  can  be  obtained.  The  location  of 
this  cutoff  point  and  the  slope  must  be  determined  for  each  new  experimental  sta¬ 
tion  by  the  methods  desr  ri bed  above .  The  total  variance  in  the  wave-height  spectrum 
can  be  determined  in  principle  from  the  coherent  component  ir.  the  received  signal  spec) 


5.5  SPECIAL  PRECAUTIONS 

In  Runs  II,  ill,  V,  VI  and  VII,  there  was  reasonably  good  agreement  between 
theory  and  experiment.  Run  IV  on  the  wther  hand  showed  o  peak  in  the  acoustic 
spectrum  not  predicted  by  theory,  whereas  Run  I  diverged  completely  from  theory. 
When  the  sou.id  spectrum  level  plots  (Chapter  4)  for  these  runs  are  examined  and 
compared  with  the  computed  curves  of  this  chapter,  it  is  seen  that  the  best  results 
are  obtained  when  the  level  of  the  sidebands  is  high  compared  with  the  level 
of  the  carrier.  For  example,  in  Run  V,  the  sidebands  are  only  about  3  dB 
below  the  carrier  and  excellent  agreement  was  obtained.  In  Run  I,  the  sidebands 
were  on  the  order  of  35  dB  below  the  carrier,  and  in  Run  IV  the  sidebands  were 
30  dB  below  the  carrier.  At  127  Hz,  the  coherent  component  dominates  the 
spectrum,  and  it  is  very  difficult  to  retrieve  useful  information.  At  this  frequency. 
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the  surface  modulation  is  very  weok  and  may  be  distorted  by  background  noise.  In 
addition,  on  very  calm  days,  such  as  the  day  Run  IV  was  made,  mixing  may  be 
reduced  in  the  water  column  and  density  gradients  that  would  invalidate  the 
assumption  of  an  isovelocity  medium  can  develop. 

It  can  therefore  be  inferred  that  the  conditions  for  obtaining  the  maximum 
information  possible  from  the  acoustic  experiment  are  as  follows: 

1.  High  acoustic  frequency  (within  the  constraints  of  the  small  roughness 
approximation). 

2.  Sufficiently  rough  surface  to  produce  vigorous  modulation  of  the  acoustic 
waves,  but  not  so  rough  as  to  violate  the  small  roughness  approximation. 

3.  Good  signal -to- noise  ratios  in  the  recording  and  processing  electronics. 

4.  The  water  column  must  be  well  mixed. 
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6.0  CONCLUSIONS 

The  surface  reradiated  signal  spectra  estimated  from  acoustic  measurements  in 
Block  Island  Sound  during  isovelocity  or  positive  gradient  conditions  are  spread  in 
frequency  and  have  two  primary  features: 

1 .  A  specular  component  exists  at  the  transmitted  frequency  and  some  scattered 
energy  exists  immediately  adjacent  to  the  specular  component.  The  scattered  energy 
effectively  broadens  the  specular  component.  This  broadening  could  be -the  result  of 
secor  d-crder  scatteringor  modulation  by  infragravity  waves,  internal  waves,  or  t*  ,es. 

2.  When  either  sea  or  swell  (but  not  both)  are  present  on  the  sea  surface, 
energy  is  contained  in  the  first-order  sidebands  of  the  acoustic  spectrum  close  to 
the  peak  frequency  in  the  surface  wave -height  spectrum.  When  sea  and  swell 
occur  simultaneously,  the  modulation  effect  of  the  swell  tends  to  mask  the  modula¬ 
tion  effects  of  the  sea. 

Second-order  sidebands  were  present  in  only  ore  measurement,  which  was 

2 

taken  at  17G2  Hz  when  the  meon-square  wave  height  was  0. 06  m  .  In  all  other 

2 

cases,  the  mean- square  wave  height  was  less  than  0. 06  m  and  no  second-order 
sidebands  co.ild  be  positively  identified. 

The  acoustic  spectra  measured  simuiJoneousiy  at  4.  5  km  and  31  km  from  the 
sound  source  are  remarkably  similar  with  regard  to  peak  location  of  firs'-order 
sidebands.  £om.  the  specular  component  and  the  sidebands  are  generally  rscrrower 
for  the  data  at  long  range.  This  occurs  because  the  high- frequency  Doppler, 
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which  corresponds  to  steep  angles  of  scatter,  is  more  heavily  attenuated  by  succes¬ 
sive  reflections  than  is  the  low-frequency  Doppler,  which  corresponds  to  small 
angles  of  scatter. 

The  scattering  is  much  stronger  at  1702  Hz  than  at  127  Hz.  For  one  experi¬ 
ment  during  nearly  similar  sea  conditions,  the  level  of  the  first-order  sidebands 
was  3  dB  below  the  carrier  level  at  1702  Hz  but  21  d&  below  the  carrier  level  at 
127  Hz. 

The  generalized  physical  optics,  long- crested  scattering  theory  predicts  that 
the  spectrum  of  an  incident  sinusoidal  signal  reradiated  by  the  rough,  moving 
surface  of  the  sea  consists  of  the  following: 

1.  A  specular  (coherent)  component  that  is  free  from  Doppler  spreading  and 
that  occurs  at  the  transmitted  (carrier)  frequency. 

2.  A  scattered  component  that  is  symmetrical  about  the  carrier;  this  compo¬ 
nent  is  the  weighted  sea  surface  height  spectrum.  The  weighting  function  depends 
on  the  signal  frequency,  the  experimental  geometry,  and  the  mean-square  wave 

height;  it  applies  only  to  conditions  of  low  surface  roughness,  i.  e. , 

-  2  ^  , 

2c  c  a  p  <  1  . 
m  n  r 

The  generalized  physical  optics,  long-crested  model  developed  in  this  study 
provides  reasonably  good  agreement  with  field  measurements,  if  the  level  of  the 
acoustic  sidebands  ere  not  more  than  25  dB  below  the  level  of  the  carrier. 

Agreement  between  theory  and  experiment  was  significantly  better  at  1702  Hz 
than  it  was  at  127  Hz.  This  follows  because  at  127  Hz,  the  effects  of  the  bottom 
are  significant  and  cannot  be  neglected.  A  normal  mode  model,  in  which  the 
surface  roughness  entered  through  the  boundary  conditions,  would  provide  better 
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results  at  this  frequency. 

The  mean-square  wave  height  can  be  determined,  according  to  theory,  from 
the  strength  of  the  coherent  component  (Eq.  42)  or  from  the  ratio  of  the  power  in 
the  sidebands  to  the  power  in  the  coherent  component.  The  ocean  wave-height 
frequency  spectrum  can  be  determined  within  a  limited  frequency  band  by  using 
measurements  of  forward-scattered  acoustic  signals  reradiated  from  the  sea  sur¬ 
face.  The  bandwidth  in  which  information  can  be  obtained  is  determined  from  the 
width  of  the  B  function,  which  depends  upon  the  geometry  of  the  experiment  and 
the  source  frequency.  More  information  can  be  obtained  from  high-frequency  sound 
than  from  low-frequency  sound  if  the  small  roughness  criterion  is  not  violated. 

This  follows  by  noting  that  the  modulation  is  wavelength  dependent  and  de¬ 
creases  with  increasing  acoustic  wavelength;  i  .e. ,  the  surface  as  seen  by  long 
acoustic  wavelengths  is  smoother  than  the  surface  as  seen  by  shor»  acoustic 
wavelengths. 

The  weighting  function  B  behaves  like  a  low-pass  filter  that  is  relatively  flat 
out  to  a  certain  point  and  then  begins  to  slope  down  to  higher  frequencies  to  a 
near-zero  or  "cutoff*'  value.  The  rate  of  slope  of  the  B  function  computed  by 
means  of  the  theoretical  model  is  greater  than  the  corresponding  slope  of  the  B 
(unction  computed  by  means  of  field  data.  This  discrepancy  results  from  the  fact 
that  the  long- crested  ocean  wave  model  does  not  contain  sufficient  energy  at 
high  frequencies,  when  the  ocean  is  really  short- crested. 

When  the  gereralized  physical  optics  theory  is  extended  to  include  the  short- 
crestedsea  surface  model  (Pierson^®),  the  slopeinfhe  B (weighting)  function  will  be 
decreased.  For  unrestricted  experimental  geomf  ry,  the  B  function  will  probably  approach 
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unity  over  the  gravity  wave  band.  Hence,  the  generalized  physical  optics  theory 
approaches  Marsh's  prediction  that  the  spectral  density  of  the  first-order  sidebands 
in  the  received,  reradiated  signal  spectrum  is  identical,  within  a  constant,  to  the 
spectral  density  of  the  ocean  wave  height. 

As  a  remote  analytical  probe  for  studying  the  dynamics  of  the  ocean  surface, 
the  acoustic  test  range  offers  exciting  possibilities,  but  much  further  study  is  still 
required.  One  feature  of  the  acoustic  method  is  that  the  acoustically  measured 
wave-height  spectrum  contains  informal  ion  about  a  very  large  area  of  the  sea 
surface  (up  to  many  thousands  of  square  meters),  whereas  the  wave  staff  is  a  point 
measurement.  Hence,  it  cannot  be  expected  that  the  spectra  estimated  from 
acoustic  and  wave-staff  methods  should  be  identical.  It  appears  that  the  acoustic 
forward-scattered  spectrum  can  provide  an  effective  means  of  studying  the  distribu¬ 
tion  of  water  wave  phase  velocities  and  the  limits  of  applicability  of  the  dispersion 
relation.  The  forward-scattered  spectrum  can  also  be  used  to  determine  the  suitability 
of  various  sea  surface  models  (such  as  the  long-crested,  the  short-crested,  and  the 
isotropic)  that  predict  the  effects  of  the  sea  surface  on  reradiated  acoustic  and 
electromagnetic  waves.  Before  such  analytical  studies  can  be  made,  however,  much 
further  research  is  indicated.  Some  topics  for  future  research  are  presented  in  the 
next  chapter. 


7.0  SUGGESTIONS  FOR  FURTHER  STUDY 


Among  the  most  difficult  problems  in  the  present  study  are  the  choice  of  size 
for  the  effective  scattering  area  on  the  ocean  surface  and  the  selection  of  me 
length  of  the  subsection  (L).  The  long-range  experimental  results  are  also  compli¬ 
cated  by  the  fact  that  Block  Island  Sound  is  a  shallow  coastal  area  with  strong  tidal 
currents;  hence,  one  is  forced  to  deal  with  data  that  have  been  contaminated  to  at 
least  some  degree  by  multipath  interference  effects. 

Since  the  surface  scattering  problem  is  rather  difficult  to  begin  with 
and  has  not  been  completely  solved  for  e^en  the  single-reflection  case,  it  does  not 
appear  wise  at  this  time  to  attempt  to  create  an  all-inclusive  theoretical  model  that 
could  predict  the  forward-scatter  spectrum  ^suiting  from  multiple  scattering  in  a 
strctified  medium  with  a  rough,  moving  boundary  and  a  boundary  with  a  variable 
acoustic  impedance.  Therefore,  with  the  main  purpose  of  this  thesis  always  in  mind 
(i.e.,  to  study  the  ocean  surface  by  means  of  reflected  and  scattered  underwater  sound 
and,  hopefully,  to  show  that  newihings  can  be  learned  ubuui  rhc  ocean  surface  spec¬ 
trum  and  the  distribution  of  phase  velocities  by  means  of  the  remote  acoustic  probe), 
a  new  experiment  is  proposed,  along  with  additional  theoretical  work. 

7. 1  THE  IDEAL  OCEAN  WAVE/ACOUSTIC 
SCATTERING  EXPERIMENT 

Although  the  omnidirectional  source  provides  the  maximum  information  about 


the  sea  surface  in  the  received  signal  spectrum,  interpretation  is  difficult.  A 
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better  means  of  obtaining  this  information,  in  which  much  greeter  control  over  the 
experiment  could  be  achieved,  would  involve  a  multi  frequency  array  and  one  or 
more  remote  receiving  arrays  in  which  the  beam  shape  and  the  angles  of  incidence 
on  the  su-face  could  be  controlled.  With  such  an  array,  the  size  of  the  scattering 
area  jould  be  easily  determined,  both  theoretically  and  by  measurements  in  the 
field.  When  the  shape  of  the  illuminated  area  is  known,  the  other  parameters  in 
l  model,  P  and  L,  could  be  precisely  determined.  In  addition,  the  array 
should  be  in  deep  water  to  eliminate  the  multipath  effects  at  low  acoustic  frequen¬ 
cies,  which  are  inevitable  in  shallow  water.  The  ocean  in  the  test  area  should  be 
as  nearly  homogeneous  as  possible  or  should  show  o  slight  increase  of  temperature 
with  depth  from  the  surface  to  the  bottom.  The  Arctic  Ocean  offers  suitable 
locations  for  the  array.  However,  if  the  theoretical  model  is  expanded  to  include 

refraction  effects,  the  area  off  Bermuda  could  be  used  and  suitable  arrays  have 
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been  installed  there.  Roderick  and  Cron  have  already  reported  on  a  surface 
scattering  experiment  conducted  in  this  area. 

With  such  a  setup,  it  would  be  possible  to  vary  the  size  of  the  illuminated 
area  on  the  surface,  to  change  the  angles  of  incidence  {with  more  than  one  receiv¬ 
ing  array),  and  to  determine  the  coherence  between  elements  of  the  receiving 
array  so  that  spatial  information  on  the  scattering  surface  could  be  obtained. 

In  addition,  much  rougher  sea  states  are  possible  ir  the  open  ocean,  and  the 
wave-height  spectrum  is  more  readily  predictable  and  is  free  from  land  mass 
refraction,  reflection,  and  diffraction  effects  that  render  basic  physical  studies 
very  difficult.  In  deep  water,  of  course,  the  dispersion  relation  is  more  nearly 


satisfied. 
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The  ocean  wave-height  spectrum  should  be  determined  by  means  of 
68 

Stillwell's  optical  technique  or  by  stereo  photography,  as  in  the  Stereo  Wave 

69 

Observation  Project  (SWOP)  experiments. 

7. 2  ADDITIONAL  THEORETICAL  STUDY 

Two  extensions  of  the  generalized  physical  optics  theory  are  evident: 

1 .  Include  the  short-crested  ocean  wave  model  (Pierson^)  to  determine 
directional  effects  and  to  compute  a  more  accurate  weighting  (B)  function. 

2.  Extend  the  theory  to  allow  for  larger  surface  roughness.  An  initial  way 
to  do  this  would  be  to  take  more  terms  in  the  power  series  expansion  for 

exp  |  2cm  c^  p  |  and  solve  the  resulting  algebraic  (but  no  longer  linear) 

equation  for  S(u)  by  numerical  means,  if  necessary.  As  a  starting  point  for  this 

generalization,  it  would  be  better  to  start  with  the  physical  optics  theory  as 
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formulated  by  Nuttall  and  Cron. 

3.  Include  the  effects  of  sound  speed  gradients  in  the  scattering  model. 


Appendix  A 

THE  ILLUMINATION  FUNCTION 


The  illumination  function  in  acoustic  scattering  theory  describes  how  the 
acoustic  intensity  is  distributed  over  the  reflecting  surface.  The  simplest  illumina¬ 
tion  function  is  the  box  car  (rectangular)  function,  which  b  1  inside  the 
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illuminated  area  and  zero  outside.  Beckmann  and  Spizz>chino  end  Perkins  °  use 

this  illumination  function. 
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Clay  and  Medwin  point  out,  however,  that  this  function  leads  to  undesirable 
diffraction  sidelobes,  which  would  be  a  serious  problem  when  acoustic  intensities 
from  many  subsections  are  summed.  To  avoid  this,  they  use  a  Gaussian  Illumination 
Function  of  the  form 

e*pj-[(x  -  XJ2/L  +  <y  -//A,] 


where 

(X  ,  Y  )  =  center  of  the  subarea 
m  m 

L  and  L^  =  lengths  of  the  subsections. 

When  the  summation  is  mode,  one  subarea  can  be  smoothed  into  the  next  by  the 
proper  choice  of  parameters. 
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in  this  study  a  more  general  method  (due  to  Nuttall  )  is  used  in  which  it  is 
possible  to  use  any  reasonable  illumination  function.  The  Gaussian  function  is 
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chossn  in  tfils  work  for  the  reasons  given  by  Clay  and  Medwin. 

The  function  H(  0,  u),  from  Eq.  (28)  in  the  muin  text,  has  been  defined  as 


/  U\  /  u  \ 

-  [  /x'n  i x+i\ 

,  u)  -  /  dx  exp(i/3x)  1 1 - 1  i  I - - J  . 

J  \  PL  /  \  PL  / 


(A-1) 


For  the  Gaussian  illumination  Function,  this  becomes 


I,  u)  =  J  dx  exp(i 


Sx)  exp  I  - 


Expanding  the  exponents  and  collecting  terms  yields 


■  wu 

'■u)  I 


2  2 

dx  exp(-2  x  +  i  0x )  exp(-  u  /2)  . 


This  integrai  can  be  evaluated  (Kinsman^)  and  the  result  is 


H(d,  u)PL{u/2)1/2  exp  j-  02P2L2/8- u2/fr2  L2  J 


(A-2) 


(A-3) 


(A-4) 


and  H  ^  d  ,  «  +  ~  / 


Now,  the  functions  H 


Eq.  (58)  in  the  main  text,  are  evaluated.  Define  the  Fourier  Transform  <5t(r)  of 
the  illumination  function  !(x)  : 


(A-5) 


A  (r)  =  — 

2r 


-co 


dx  exp(-irx)  Ifc)  ; 
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therefore, 


+oo 


I  to 


-  J  dr  exp(i 


p(irv)  A  (r)  . 


(A-6) 


-oo 


Substituting  Eq.  (A-6)  into  Eq,  (54)  of  the  main  text  yields 


+oo 


-00 


H  a  -  =  J  ^ °  6X^  ”*  U  J  dxexp(i0x) 


+00 


-CD 


*  jj  do'  db'  expj  i  ^  (“  )  ] 

-00  ( 


a(q)  A(b) 


Factoring  the  terms  in  Eq.  (A-7),  interchanging  the  order  of  integration,  and 
collecting  terms  in  u  and  x  result  in 

+oo 

‘  |e*p 


(A-7) 


/  djexTi[(t''7)+°^b]u 


-CO 


+oo 


+ao 


+oo 


J  dx  exp  i  +  ° — ^  J  x  |  da  &  (a)  J  d\>  $  (l>) 


-co 


(A-6) 


-oo 


-ao 
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The  result  of  integrating  Eq.  (A-8)  with  respect  to  u  and  x  is 

+oo  +tjo 

a “  J  da'|^(a)  J  dl>  £ (i>)  2*S  j^a  «  ^-^+  ° — »"j 


-00  -CD 


•2*5 


&  + 


/ .  ./'ll 
a  +  b  5 1 

PL  Jj' 


£4-9) 


Collecting  terms  and  integrating  with  respect  to  b  yield 

+0D 

H  ^/3,  a  ~  ~~  f  d°  ^  ^  P +  ^  (a  "  j  ^ 


2  2  2 
PL 


-ao 


•  25|2a  +  PL/3  +  2PL 


m- 


(A— 10) 


Integrating  Eq.  (A-10)  with  respect  to  a  yields 


h(p,  4hP2L2  £  PL/3/2 -PL^a-^J 


(A-ll) 


For  the  Gaussian  Illumination  Function  the  result  is 
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+oo 

J-PL0/2  -  PL  fa-  sJ \]=”  f  dx  exp-i|-PL(a-— \  -PL  0/2 

L  \  g/j  2ii  (  L  \  g/ 

-oo 


X  -  X 


(A-12) 


+00 


[jJ.pL0/2  +  PL^a-l-J  =—  J  dx  exp  -ij^PL  ^a-^-PL^jx  -  x2| 


-00 


When  £q.  (A-12)  is  Integral  using  the  method  of  Kinsman., 


the  result  is 


S  -PL0/2  -PL 


2» 


w 


+  p2lV/4 


(A- 13) 


Substituting  Eq.  (A-13)  into  Eq.  (A-ll)  yields  the  final  result  for  the  function 


in  the  case  of  the  Gaussian  Illumination  Function: 
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■K) 


B2,2 

P  L  exp 


A  similar  procedure  is  followed  to  obtain 


in  h|/>,  <*  +  — j 


(A-14) 
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Appendix  B 

OCEANOGRAPHIC  OBSERVATIONS 
DURING  THE  BLOCK  ISLAND  SOUND  TESTS 

Hourly  weather  observations  for  the  time  period  covering  the  acoustic/oceano- 
graphic  experiments  in  Block  Island  Sound  were  furnished  by  William  B.  Phelan,  of 
the  U.  S.  Weather  Bureau.  The  data  were  taken  at  Block  Island  Airport  and  are 
presented  in  Figure  B-1. 

Sound  speed  versus  depth  data  obtained  during  the  measurements  are  given  in 

Table  B-1.  The  values  were  computed  from  temperature  and  salinity  data  by  means 
25 

of  Wilson's  equation. 


TABLE  B-l 


SOUND  SPEED  VERSUS  DEPTH 


Date  Site 

•  ime 
(hr) 

Depth 

(ft) 

Sound  Speed 
(m/sec) 

29  December  1969  Station  B 

1025 

5 

1470. 10 

40 

1469.90 

80 

1470. 10 

114 

1469.80 

Station  C 

115t> 

5 

1469.20 

40 

1469.00 

80 

1469.20 

113 

1469.40 

Station  A 

1717 

5 

1466. 15 

54 

1466.40 

107 

1467.40 

160 

1467.65 

8  January  1970  Station  B 

1005 

5 

1464.75 

40 

1465.35 

80 

1466.05 

114 

1466.25 

Station  C 

1107 

5 

1465.35 

35 

1465.25 

72 

1463.60 

110 

1464.20 

Station  A 

1640 

5 

1459.80 

35 

1460.00 

75 

1459. 70 

110 

1459.50 

From  Nalwalk,  Rath  bun,  Robinson, 

and  Riley.  ^ 

TABLE  B-l  (Confd) 
SOUND  SPEED  VERSUS  DEPTH 
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Date 

Site 

Time 

Depth 

Sound  Speed 

(hr) 

(ft) 

(m/sec) 

16  January  1970 

Station  B 

0945 

5 

1455. 25 

39 

1455.45 

78 

1456. 15 

113 

1459.45 

Station  C 

1035 

5 

1455. 75 

40 

1456. 75 

80 

1457.65 

114 

1460.65 

Station  A 

1515 

5 

1456. 15 

43 

1456.65 

86 

1456.85 

124 

1457.25 

26  January  1970 

Station  B 

1000 

5 

1452.20 

40 

1454.65 

80 

1457.35 

120 

1457.65 

Station  C 

1102 

5 

1453.80 

37 

1455.65 

75 

1456.65 

115 

1456.85 

Station  A 

1710 

5 

1451.20 

33 

1452.60 

66 

1452.40 

97 

1452.80 
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TABLE  B-l  (Cont'd) 

SOUND  SPEED  VERSUS  DEPTH 

Date 

Site 

Time 

(hr) 

Depth 

(ft) 

Sound  Speed 
(m/sec) 

30  January  *970 

Station  B 

0955 

5 

1452.60 

40 

1453.60 

80 

1449.85 

116 

1457.60  i 

j 

Station  C 

1110 

5 

1453.50 

39 

1453.80 

78 

1455.20 

113 

1455.40 

Station  A 

1346 

5 

1453.30 

34 

1453.50 

68 

1454. 10 

97 

1454.80 

- i 

Figure  B-1.  (Cont*d)  Hourly  weother  observations  at  the  Block  Island  Airport 


Figure  B-l.  (Cont'd)  Hourly  weather  observations  at  the  Block  Island  Airport 
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Figure  B-l.  (Cont*d)  Hourly  weather  observations  at  the  Block  Island  Airport 
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APPENDIX  C 

INSTRUMENTATION 

The  following  major  instrumentation  was  used  in  this  study: 

1.  An  XU-1378  magnetostrictive  scroll  transducer.  —The  XU-1378  has  c 
resonance  frequency  of  1700  Hz,  a  bandwidth  of  200  Hz,  and  a  source  level  of 
+97.6  dB//l  |ibar.  It  is  connected  to  shore  via  cables  that  terminate  at  the  Block 
Island  field  station. 

2.  A  127- Hz  Minneapolis- Honeywell  source.  —This  projector  has  a  band¬ 
width  of  25  Hz  and  a  source  level  of  +106. 2  dB//l  pbar  and  consists  of  one  stave 
of  an  HX-37  ceramic  bender-bar  projector.  During  these  experiments,  one  of  the 
three  stages  of  the  projector  failed;  consequently,  the  source  level  was  somewhat 
less  than  the  output  specified  for  three  stages. 

3.  Two  bottom- mounted  DT-55  hydrophones.  —The  hydrophones  were 
implanted  about  1-m  above  the  bottom  south  of  Fishers  island  and  were  connected 
by  cable  to  shore.  Each  of  the  hydrophones  has  a  self-contained  preamplifier. 

4.  An  NUS  LM- 3  hydrophone.— The  LM-3  contains  preamplifier  and 
calibration  circuitry.  It  is  omnidirectional  and  completely  flat  in  response  from 
10  Hz  to  100  kHz.  The  sensitivity  is  -90.4  dBV//l  pbar.  This  hydrophone  was 
connected  by  cable  to  the  ship's  laboratory. 

5.  A  standard  Navy  model  Bendix  anemometer.  —This  device  was  mounted 
on  top  of  the  Block  island  field  station. 


6.  A  Broincon  Type  474  wove  measuring  system.  —  This  system  consists  of  a 
free-floating  accelerometer  wave  buoy.  Model  231-2,  tethered  to  a  toroid  moor¬ 
ing  buoy  containing  a  transmitter  and  batteries.  The  buoy  can  be  turned  on  from 
shore.  The  data  are  received  at  the  Fishers  Island  field  station.  The  shape  of  the 
buoy  is  c  truncated  cone  with  a  small  cylinder  extending  down  from  its  base.  In 
the  normal  position,  this  configuration  provides  a  flat  undersurface  with  good  wave- 
following  characteristics.  The  buoy  is  made  of  fiberglass- reinforced  plastic  and  is 
self-righting. 

The  principle  employed  in  the  accelerometer  is  that  of  measuring  the  rise  and 
fall  time  of  a  moving  body  (a  coil),  which  is  given  an  upward  impulse  of  fixed 
momentum  relative  to  a  reference  frame.  The  coil  is  located  in  an  annular  air  gap 
in  a  permanent  magnet.  When  contacts  close  at  the  base  of  the  gap,  a  monostable 
muitivibratoi  circuit  supplies  a  current  pulse  of  fixed  amplitude  and  duration.  The 
system  is  essentially  pulse-code  modulated,  i.e.,  the  varying  acceleration  due  to 
the  ccean  waves  modulates  the  basic  pulse  repetition  rate. 

7.  A  12-m  aluminum  spar  buoy.  —The  buoy  was  designed  and  built  by  Kirk 
Patton,  of  the  Naval  Underwater  Systems  Center  (NUSC). 

8.  An  Interstate  Electronics  Model  438-2000  wave  staff.— The  staff  was 
constructed  from  1/2-in.  PVC  pipe  and  28-gauge  tophat-A  resistance  wire, 
which  served  as  the  active  element.  The  connection  to  the  end  of  the  resistance 
w;re,  which  is  normally  in  the  sea  water,  is  a  length  of  20-gauge,  vinyl- insulated 
wire.  This  insulated  wire  is  recessed  into  a  longitudinal  slot  in  the  PVC  pipe  and 
is  held  in  place  by  the  resistance  wire,  which  is  spirally  wound  in  grooves  in  the 
PVC.  The  15- ft  wave  staff  has  a  standard  resistance  of  44  ohms/linear  foot  and  is 
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clamped  at  both  ends  to  the  aluminum  spar  buoy  so  that  the  staff  proper  is  6  in. 
away  from  the  buoy. 

The  staff  was  powered  by  a  stable,  dc  power  supply  at  12  V  and  20  mA.  The 
staff  output  (0-5  Vdc)  is  proportional  to  wave  height.  The  calibration  accuracy  is 
1%  of  the  dynamic  range,  and  the  linearity  is  1%  over  90% of  the  dynamic  range. 

The  voltage  to  the  sensor  unit  is  5  VRMS  at  500  Hz.  The  electronics  package 
contains  fully  transistorized  circuitry  and  is  housed  in  a  PVC  cylinder.  A  cali¬ 
bration  check  performed  at  the  NIISC  Millstone  Pond  facility,  along  with  a  block 
diagram  of  the  system,  is  shown  in  Figure  C-l.  A  system  calibration  was  performed 
in  the  field  before  and  after  each  run. 

9.  An  NUS  TR-4  velocimeter.  —The  NUS  TR-4  consists  f  a  pair  of  piezo¬ 
electric  ceramic  transducers  and  an  acoustic  reflector  rigidly  mounted  to  form  a 
sound  path  of  known  geometry.  The  sound  path  is  folded  to  minimize  Doppler  error 
due  to  fluid  flow  along  the  acoustic  path.  The  on-deck  equipment  consists  of  an 
amplifier,  a  frequency  doubler,  and  an  electronic  counter.  The  frequency  of  the 
counter  was  read  at  3-m  depth  intervals  and  was  converted  to  sound  speed  by  an 
equation  of  the  form  C  =  Kf,  where  C  is  the  sound  speed,  f  is  the  frequency, 
and  K  is  the  constant.  The  accuracy  of  this  instrument  is  0.  15  m/sec  with  a 
stability  of  0. 01  m/sec. 

10.  General-Radio,  coherent  decade  frequency  synthesizers  (Type  1161). 

—  One  synthesizer  was  used  to  control  the  output  frequency  of  the  sound  source. 

The  other  synthesizers  were  used  to  generate  the  reference  signals  that  were  recorded 
on  the  magnetic  tape  with  the  data  signal  to  provide  for  tape  speed  control  and 
band  shifting.  This  instrument  is  a  tunable  quartz-crystal -control  led  oscillator 
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with  resolution  to  nine  significant  figures.  It  has  an  rms  fractional  frequency 

9 

deviation  on  the  order  of  3  parts  in  10  . 

Although  these  instruments  provide  an  extremely  accurate  and  stable  output,  it 
was  not  always  possible  to  set  the  frequency  difference  between  the  transmitted 
signal  and  the  center  frequency  of  the  received  signal  to  exactly  8  Hz.  This 
difficulty  occurred  because,  on  some  occasions,  frequency  counters  with  sufficient 
accuracy  to  check  the  output  of  the  synthesizers  were  not  available.  Hence, 
although  no  drift  occurred  in  the  synthesizer  output,  the  absolute  frequency  was 
slightly  in  error.  A  central  timing  control  system  could  have  been  used  to  lock  all 
synthesizers  at  the  desired  frequency,  but  none  was  available.  As  indicated  in  the 
text,  this  error  is  corrected  by  measuring  frequency  shifts  in  the  received  signal 
from  the  peak  observed  carrier  frequency,  rather  than  from  8  Hz. 

1 1.  A  Montronics  frequency  synthesizer.  —  This  synthesizer  was  used  to  record 
a  stable  frequency  on  the  magnetic  tape  for  use  in  analog- to- digital  conversion  of 
the  R/V  UCONN  data. 

12.  An  Ampex  Model  FR-1200,  14-trock,  analog  tape  recorder.— The  signals 
from  the  Fishers  Island  hydrophone  were  FM  recorded  at  7  1/2  IPS.  At  this  tape 
speed,  the  frequency  response  to  within  1  dB  is  0  to  2500  Hz,  with  a  signa'i-to- 
noise  ratio  of  42  dB  (broad  band)  and  a  total  harmonic  distortion  of  1. 5%.  The 
record- reproduce  voltage  linearity  is  1.0%  of  full  band. 

13.  An  Ampex  Model  CP  10QD,  7  track,  onalog  tope  recorder  on  R/V 
UCONN.  —  The  data  were  FM  recorded  at  7  1/2  IPS.  At  this  tape  speed,  the 
frequency  response  to  within  1  dB  is  0  to  2500  Hz,  with  a  broad- band  signal -to- 
noise  ratio  of  42  dB  and  a  1-Hz  band  signal -to- noise  ratio  of  60  dB.  Tape  flutter 
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is  rated  at  less  than  0. 25%. 

14.  Two  Intronics  Model  M502  wideband  analog  multipliers.— The  multipliers 
were  employed  to  bandshift  the  acoustic  and  radio  signals  to  8  Hz  prior  to  recording. 
These  instruments  feature  four  quadrant  dc  operation  on  both  inputs  and  employ 
variable  transconductance  to  maintain  a  smooth  linearity.  The  full-scale  accuracy 
is  +0.  5%  in  amplitude.  Maximum  deviation  from  linearity  is  + .  3%. 

The  instrumentation  used  during  the  accelerometer  buoy  experiments  is  shown 
in  Figure  C-2.  The  acoustic  signal  was  recorded  on  tape  and  later  bandshifted  to 
2  Hz  for  computer  processing.  The  accelerometer  buoy  signal  was  played  into  a 
firequency-to-voltage  converter,  the  output  of  which  was  digitized  for  computer 
analysis. 

For  Runs  1  through  7,  the  acoustic  signal  was  bandshifted  to  S  Hz  before  it  was 

recorded  on  magnetic  tape.  Extreme  care  was  taken  to  eliminate  the  effects  of 

tape  flutter.  A  12. 5  kHz  signal  generated  by  a  quartz-crystal,  frequency  synthe- 

71 

sizer  (Montronics)  was  recorded  and  later  used  in  a  servocontrol  phase  lock 
system  on  playback  during  analog- to -digital  conversion.  In  addition,  a  256-Hz 
sampling  frequency  from  a  stable  oscillator,  as  well  as  voice  and  time  code,  was 
recorded.  Block  diagrams  of  the  electronics  aboard  R/V  UCONN  and  at  the 
NUSC  laboratory  are  shown  in  Figures  C-3  and  C-4,  respectively. 

In  a  later  version  of  these  experiments,  performed  in  late  January  1970,  the 
acoustic  signal  was  transmitted  over  a  radio  link  and  the  bandshifted  radio  frequency 
was  ^corded  simultaneously  with  the  bandshifted  acoustic  signal  and  the  wave  staff 
output,  as  shown  in  Figure  C-5.  The  purpose  of  this  modification  was  to  prove 
that  none  of  the  sideband  structure  seen  around  the  acoustic  signal  spectrum  was 
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due  to  instrumentation.  Some  other  minor  changes  were  also  made  to  improve  the 
quality  of  the  recording,  such  as  using  a  bucking  voltage  to  eliminate  the  dc  com¬ 
ponent  of  the  wave  staff.  In  addition,  the  internal  13. 5- kHz  FM  subcarrier  was 
recorded  and  used  on  playback  in  a  new  flutter  compensation  system,  which 
improved  the  signal -to- noise  ratio  by  6  dB.  ^ 
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Figure  C-2.  Recording  electronics  for  the  accelerometer  buoy  experiment 
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Figure  C-3.  Receiving  and  recording  electronics  aboard  R/V  UCONN 

(29  December  1969  -  25  January  1970) 
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Figure  C-4.  Recording  electronics  at  the  NavaF  Underwater  Systems  Center 

(29  December  -  30  January  1970) 


WAVE  STAFF 
(0-5Vdc) 
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Figure  C-5.  Modified  receiving  and  recording  electronics  aboard  R/V  UCONN 

(26  January  -  30  January  1970) 
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Appendix  D 

DATA  PROCESSING  AND  ANALYSIS 

ANALOG  PROCESSING  AND  ANALOG- 
TO -DIGITAL  CONVERSION 

Data  taken  in  the  field  were  FM  recorded  on  analog  magnetic  tape,  and 
Krohn-H'Ie  low-pass  and  band-pass  filters  were  used  to  suppress  the  high-frequency 
content  out  of  the  band  of  interest.  On  playback  during  analog- to-digital  (A-D) 
conversion,  additional  low- pass  filtering  prevented  future  aliasing  of  the  spectra. 

A  Control  Data  Corporation  (CDC)  A-D  converter  digitized  the  wave-staff  and  the 
bandshifted  acoustic  data  at  a  rate  of  64  samples  per  second,  whereas  the  ambient  noise 
data  and  calibration  signals  were  digitized  at  a  rate  of  8192  samples  per  second. 

The  control  signals  from  the  frequency  synthesizers,  the  256-Hz  sampling 
frequency,  and  the  FM  subcarrier  shorted  input  were  recorded  on  magnetic  tape; 
during  playback,  these  signals  were  used  to  eliminate  the  effects  of  tope  flutter 
from  the  recorded  data  (P.  Stahl Digital  data  were  recorded  on  digital  com¬ 
puter  tapes  in  CDC  format  for  further  computer  processing. 

METHODS  OF  ANALYSIS 

In  studies  involving  spectrum  analysis  of  digital  data,  resolution  and  stability 
are  primary  considerations.  Resolution  determines  the  ability  to  observe  fine 
detail  in  the  spectrum,  whemas  stability  is  a  measure  of  the  variance  of  the 
smoothed  spectral  estimator.  Resolution  and  stability  determine  the  length  of  the 


data  record  required  to  answer,  with  confidence,  the  questions  being  asked  about 
the  spectrum. 

In  general,  it  is  always  an  advantage  to  obtain  as  long  a  data  record  as  possible. 

However,  the  process  being  observed  must  be  reasonably  stationary  over  the  length 

of  time  of  the  observations.  For  wind  waves  in  Block  Island  Sound,  this  time  is  on 

the  order  of  20  min,  which  provides  an  upper  bound  on  the  length  of  data  recorded 

for  subsequent  processing  to  determine  spectral  content. 

In  addition  to  considerations  of  resolution  and  stability,  Bingham,  Godfrey, 

58  57 

and  Tukey  and  Welch  point  cut  that,  prior  to  spectrum  analysis,  the  linear 

trend  and  dc  component  of  the  data  should  be  removed.  Spectral  leakage  must  be 

minimized  prior  to  scjuorinQ  end  adding  the  Fourior  co€?ficisnts.  Ai*  oppropriot© 

spectral  window  must  be  chosen  for  this  purpose.  Statistical  stability  in  the 

estimates  computed  by  means  of  the  Fast  Fourier  Transform  (FFT)  can  be  obtained 

57 

only  by  averaging  several  estimates  of  the  spectrum.  Welch  has  developed  a 

method  of  averaging  that  employs  overlapped  sections  of  data  to  achieve  maximum 

variance  reduction  for  a  given  resolution  and  record  length. 

72 

These  considerations  have  been  utilized  by  Nuttall  to  develop  an  optimum 
FFT  spectral  analysis  procedure,  which  has  been  used  to  process  the  data  presented 
in  Chapter  4. 

in  order  to  obtain  the  data  presented  in  Chapter  4,  the  following  essential 
steps  were  followed: 

1.  Determine  the  resolution  desired.  This  specifies  the  length  of  each  piece 
of  record  to  be  used  in  computing  the  average  spectrum. 

2.  Remove  the  dc  component  and  the  linear  trend  from  each  piece. 


3.  Window  the  data  in  the  time  domain.  The  window  used  on  Runs  l-VN  is 
the  cosine  window,  which  reduces  the  maximum  sidelobe  level  to  -31.5  dB.  The 
recommended  overlap  for  the  cosine  window  is  60%.  However,  computational 
convenience  dictated  an  overlap  of  50%.  For  the  cosine  window,  the  resolution 

^3dB  ^  ^  's^ 


B^g  =  1.44/segment  length 


(D-l) 


4.  Compute  the  FFT  on  each  piece.  Square  and  add  to  determine  the 
"power"  and  then  form  the  average.  The  result  is  a  stable  estimate  of  the  spectrum 
if  a  sufficient  number  of  pieces  have  been  employed  in  the  average. 

The  equivalent  degrees  of  freedom  (EDF)  is  defined^  as 


EDF  = 


Var 


M ' 


(P-2) 


where 


E  |Sp(f) ^  =  expected  value  of  the  estimate  of  the  true  spectrum  S(f)  at 
frequency  f ,  based  on  p  pieces. 

Var  ”  variance. 
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For  the  cosine  window  for  50%  overlap,  Nuttall  has  shown  that 


EDF  * 2-  76«W  - 


(0-3) 


where  T  =  total  record  length.  The  number  of  pieces  into  which  the  record  must 


be  sectioned  to  obtain  this  number  of  EDF  is 
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EDF  =  >-39  8^1-1  .  (CM) 

The  data  discussed  in  Chapter  4  were  spectrally  analyzed  by  means  of  NUSC 
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Computer  Program  S 1086/W  (Arnold,  Nuttall,  Ferrie,  and  Carter  )  in  order  to 
accomplish  the  above  procedures.  Several  resolutions  were  used  in  the  analysis, 
as  suggested  by  Jenkins  and  Watts.  ^  The  resolutions  were  0. 0056  Hz,  0. 01 12  Hz, 
0. 0225  Hz,  and  0. 0450  Hz.  These  resolutions  correspond  to  transform  as  of 
8192,  4096,  2048,  and  1024  data  points,  respectively  The  90%  confidence 
intervals  were  computed  by  means  of  the  approximating  chi-square  distribution,  as 
given  by  Blackman.  /'i  The  90%  confidence  limits  are  approximately  equal  to 

10  [ - ! -  +  - L_ |  dB  for  EDF  >3  .  (D-5) 

\  3  EDF  -1  “  (EDF  -  1)  '/ 

AN  EXAMPLE  OF  A  SPECTRUM  ANALYSIS 

In  the  main  text,  the  acoustic  signal  record  for  Run  V,  at  Position  C  was 
nearly  20  min  long.  The  sampling  rate  during  A-D  conversion  was  64  Hz.  Only 
every  other  data  point  was  used  in  the  spectrum  analysis  in  order  to  permit  finer 
resolution  than  would  be  possible  if  every  data  point  were  utilized.  This  results  in 
a  time  increment  between  data  samples  of  0. 03125  sec. 

Since  high  resolution  was  desired,  the  maximum  transform  size  permitted  by  the 
Cooley-Tukey  algorithm^  on  the  UNIVAC  1 108  was  tried  first.  Figure  D-l  shows 
the  result  of  the  calculations.  The  -3  dB  equivalent  resolution  bandwidth  for  this 
case  is  0.0056  Hz.  However,  with  only  10  degrees  of  freedom,  the  90  7o  confidence 
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Intervals  an  very  large  and  thus  result  in  poor  statistical  stability.  Clearly,  this 
record  has  been  over -resolved. 

Consequent*}',  the  transform  size  was  halved,  and  the  spectrum  analysis  was 
repeated  (Figure  D-2).  Almost  twice  as  many  pieces  were  averaged,  and  the 
number  of  degrees  of  freedom  increased  to  21,  However,  the  bandwidth  increased 
to  0.0112  Hz.  The  confidence  interval  was  still  too  large,  and  it  appears  that  this 
speotruut  was  a  I  sc  over-resolved. 

Figure  D-3  shows  the  results  on  the  same  record  calculated  for  a  transform  size 
of  2048  data  points,  which  permits  twice  as  much  averaging  as  in  the  previous  ease. 

The  bandwidth  is  still  satisfactory  (0.0225  Hz),  whereas  the  confidence  interval  is 
reasonably  small.  This  estimate  of  the  spectrum  has  45  degrees  of  freedom. 

The  analysis  was  carried  one  step  further,  i.  e. ,  a  transform  size  of  1024  points 
was  used  to  generate  an  average  spectrum  with  91  degrees  of  freedom  (Figure  D-4). 

This  estimate  of  the  spectrum  has  excellent  stability  but  does  not  reveal  much  de¬ 
tail.  The  bandwidth  is  0.0450  Hz,  which  is  too  coarse  since  ocean  swell  peaks  can 
occur  at  0.08  Hz.  From  this  andother  examples,  it  is  concluded  that  a  resolution  of 
0.0225  Hz  (a  transform  size  of  2048  points)  provides  the  optimum  estimate  of  the 
spectrum  for  the  purposes  of  this  study. 

Program  51086/T  was  written  by  J.  F.  Ferrie,  of  NUSC,  to  compute  the 
acoustic  and  wave-height  spectra  by  means  of  S1086/W  and  output  the  result  on 
magnetic  tape.  These  tapes  were  used  as  input  to  the  theoretical  Program  SI 731. 

Since  the  wave  staff  is  not  precisely  linear  (Figure  D-1),  the  spline  interpola¬ 
tion  routine  of  E.  Mehr  (New  York  University,  Department  of  Meteorology  and 
Oceanography)  was  used  to  convert  from  wave-svaff  voltage  to  wave  height. 


H ir*&r*a  *W3JHJ^/S5WI  ff^ffp*! 


Figure  D-4.  Spectrum  level  of  ocoustic  signal  at  Position  C 
(resolution  =  0. 0550  Hz,  EDF  =  91) 
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COMPUTER  PROGRAM  FOR  CALCULATING  THE 
SPECTRA  FROM  THE  PHYSICAL  OPTICS  MODEL 


Computer  Program  SI  731,  which  is  included  in  tho  foiiowffljji  pages,  was  written 
by  R.  C.  Jennings,  of  NUSC,  to  calculate  the  acoustic  arwlsbcean  wave-height 


n^ce 


spectra  by  means  of  Eqs.  (67)  and  (69),  respectively.  The  program  is  written  in  the 
FORTRAN  V  language  for  the  UNIVAC  1108  Computer  and  the  Stromberg-Cos  lson 
406P  Plotter.  The  input  to  this  program  is  the  tape  output  from  the  Spectral  Analysis 
Program  S  108£/l  (Appendix  D). 

The  output  of  Program  S1731  includes  the  following  computer-generated  plots: 

1.  The  computed  weighting  function  B(f)  versus  f. 

2.  Measured  B(f)  versus  f. 

3.  Aj  (f)  measured  and  A ^ (f)  calculated  versus  f. 

4.  S(f)  measured  and  S(f)  calculated  versus  f. 

Figures  33  through  43,  in  Chapter  5  of  the  main  text,  are  outputs  of  this  program. 
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